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>• ' Abstract 

o: 

We present a detailed analysis of the structure of the conservation laws in quantum 
f — ■ integrable chains of the XYZ-type and in the Hubbard model. The essential tool for the 

former class of models is the boost operator, which provides a recursive way of calculation of 
the integrals of motion. With its help, we establish the general form of the XYZ conserved 
charges in terms of simple polynomials in spin variables and derive recursion relations for 
the relative coefficients of these polynomials. Although these relations are difficult to solve 
tj- \ in general, a subset of the coefficients can be determined. Moreover, for two submodels of 

the XYZ chain - namely the XXX and XY cases, all the charges can be calculated in closed 
form. Using this approach, we rederive the known expressions for the XY charges in a novel 
way. For the XXX case, a simple description of conserved charges is found in terms of a 
Catalan tree. This construction is generalized for the su(M) invariant integrable chain. We 
also investigate the circumstances permitting the existence of a recursive (ladder) operator 
in general quantum integrable systems. We indicate that a quantum ladder operator can 
be traced back to the presence of a hamiltonian mastersymmetry of degree one in the 
classical continuous version of the model. In this way, quantum chains endowed with a 
recursive structure can be identified from the properties of their classical relatives. We also 
show that in the quantum continuous limits of the XYZ model, the ladder property of the 
boost operator disappears. For the Hubbard model we demonstrate the non-existence of a 
ladder operator. Nevertheless, the general structure of the conserved charges is indicated, 
and the expression for the terms linear in the model's free parameter for all charges is 
derived in closed form. 

* Work supported by NSERC (Canada). 
10/94 
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1. Introduction 

This work is concerned with a detailed analysis of the structure of the conservation 
laws for quantum integrable chains with short range interactions. Attention is focused on 
the XYZ model: 

H = ^2[\ x <Tj tf+1 + A y (rJ(jJ +1 + \ x <Tj(Tj+l\, (1-1) 

jeA 

(where A is the spin lattice, a^ ,y ' z are the Pauli sigma matrices, acting non-trivially only 
on the i-th site, and X x , X y , X z are constants), the Hubbard model: 

H = EW°i+i + "K+i + *W+i + T W + i + U ^ C 1 - 2 ) 

(with cr and r standing for two independent sets of sigma matrices, and U a constant), 
and their close relatives. A brief survey of known results precedes the presentation of our 
findings. 

1.1. The problem of calculating conservation laws and ladder operators in integrable sys- 
tems 

Beyond the existence proofs, very few explicit results are known on the structure of 
the conserved charges of integrable quantum chains, even for rather basic systems such as 
the XYZ or the Hubbard model. Most of the results concerning these two particular classes 
of models were obtained by means of the Bethe ansatz method, whose immediate concern 
is the construction of the hamiltonian eigenstates (see [1] for the XYZ-type models, and 
[2] for the Hubbard model). For the XYZ model, a second input came from the discovery 
of its equivalence with an exactly solvable two-dimensional classical statistical model, the 
eight- vertex model [3]. The row-to-row transfer matrix of the latter model, denoted by 
T(A), where A is a spectral parameter arising in the parameterization of the Boltzmann 
weights, turns out to commute with the XYZ hamiltonian [4]: 

[T(A),ff]=0. (1.3) 

The expansion of the logarithm of the transfer matrix in terms of the spectral parameter 
yields a family of conserved charges. This family is infinite in the infinite chain limit. 
Furthermore, the row-to-row transfer matrices evaluated for different values of the spectral 
parameter were shown to commute [5]: 

[T(A),2ty)] = 0. (1.4) 
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This directly implies the mutual commutativity of all these conserved charges and thereby 
proves the complete integrability of the XYZ chain. 

The last statement embodies implicitly the standard working definition of integrability 
for a quantum system with N degrees of freedom, namely the existence of N independent, 
mutually commuting integrals of motion. This aspect (the existence of a family of commut- 
ing charges) of a system that can be solved by the Bethe ansatz in its original (coordinate) 
version is not immediate, although non-questionable. In the context of the Hubbard model, 
the absence of a direct link between the coordinate Bethe ansatz and integrability, moti- 
vated Shastry [6] to parallel the integrability proof of the XYZ model by constructing the 
two-dimensional classical statistical version of the Hubbard model, whose transfer matrix 
is related to the Hubbard hamiltonian. The demonstration of (1.4) for this transfer matrix 
establishes then "directly" the integrability of the Hubbard model. It should be added 
that in the quantum inverse scattering reformulation of the Bethe ansatz (the algebraic 
Bethe ansatz) developed in the last 15 years [7, 8], the transfer matrix is a central object, 
which makes algebraic Bethe ansatz solvability and quantum integrability in the above 
sense, essentially equivalent. However, the Hubbard model has not yet be reformulated 
from that point of view. 

Recent works on non-abelian symmetries in quantum chains call for a clarifying com- 
ment concerning the relation between the transfer matrix and the conserved charges. The 
expansion of the logarithm of the transfer matrix in terms of the spectral parameter at 
A = (or some other suitable finite value) leads to local conserved integrals of motion, 
which all mutually commute [9] . The first member of this infinite sequence is proportional 
to the defining hamiltonian of the quantum chain. Locality means that interaction involv- 
ing a certain set of sites disappears when the distances between them become sufficiently 
large. On the other hand, the expansion of the monodromy matrix at A = oo produces 
conserved integrals which are non-local and non- commuting (see e.g. the review [10] which 
focuses on the isotropic case, i.e. the XXX model). More precisely, these charges all 
commute with the hamiltonian and the higher local charges, but not among themselves. 
They in fact generate a highly non-trivial non-abelian algebra. For the XXX model for 
instance, this algebra provides a realization of a Yangian [11]. Notice moreover that the 
commutativity of these non-local charges with the local ones is true only in the limit of 
an infinite chain, in contrast with the local ones, which are conserved also for finite chains 
with periodic boundary conditions. In this work, we are mainly interested in the local 
charges and we will have little to say about the non-local ones. 



The most direct way of proving the integrability of a quantum chain, independently 
of (any) Bethe-ansatz formalism or a relation with a two-dimensional statistical system, is 
certainly to display a quantum Lax equation that reproduces the hamiltonian equation of 
motion of the chain. This construction can be done rather systematically. For the XYZ 
chain and the Hubbard model, the results can be found respectively in [12, 13] and [14]. 
However, this framework does not lead to a new way of calculating the conserved integrals. 
Indeed, the Lax operator depends upon a spectral parameter and the conserved laws are 
recovered from a series expansion, in terms of this spectral parameter, of the product of 
all the Lax operators, one for each site. This product is nothing but the transfer matrix. 

Therefore, at our present state of knowledge, a systematic computation of the con- 
served integrals in a generic quantum integrable chain must necessarily proceed through 
the series expansion of the generating function, i.e. the transfer matrix. But for actual 
computations, this method is completely impractical, apart from the evaluation of the 
hamiltonian itself and maybe a first few non-trivial conserved laws. The reason for this is 
that the transfer matrix is a formidable object, with its size growing exponentially with 
the number of spins. Even relatively small chains (|A| ~ 10) provide a difficult challenge 
for computer algebra programs, requiring vast amounts of running time and memory. 

A different approach is thus clearly desirable. Fortunately, in the XYZ model there 
exist an alternative, recursive way of calculating the conserved charges. It is based on the 
existence of a ladder operator, i.e. an operator B with the property: 

[5,Q n ]~Q n+ i, (1.5) 

where Q n denotes a charge with at most n adjacent spins interacting, with Q2 being 
proportional to the hamiltonianJil For the XYZ model, such an operator is provided by 
the first moment of the hamiltonian density (see e.g. [15, 16]); i.e., if one writes the 
hamiltonian in the form 

H ~ £>iJ+i> (1-6) 

jeA 

it reads 

* = Z>* *ij+i- (1-7) 



Note that we allow for the possibility of a linear combination of lower order charges Q m <r 
on the (rhs) of 1.5. 



The operator B generates Lorentz transformations, by acting on the transfer matrix 
by differentiation with respect to the spectral parameter [17, 18]: 

[*,T(A)] = ^l (1.8) 

Due to (1.8), B is usually called a boost operator. For models with short-range interaction, 
which are related to a statistical system via the transfer matrix formalism, both (1.8) and 
the ladder property (1.5) are equivalent. But since we emphasize the latter property, we 
will most often use the more appropriate ladder qualitative. 

The formula (1.5) is a useful tool. Not only does it provide a convenient recursive 
way to calculate explicit expressions for the conserved charges, avoiding the computational 
difficulties with the expansion of the transfer matrix, but, as we will see below, it is also 
helpful in proving general properties of the family {Q n }. However, even with the use of 
(1.5), the complexity of calculations of higher order charges is still considerable. Moreover, 
the problem is not only in the computational complexity, but also that the results, at first 
sight, seem hopelessly complicated and no pattern can be easily discerned. In one of the 
earlier works concerned with calculations of the explicit forma of the conservation laws 
[21], the charge Qq is described in the following way: [it] "is really a monster representing 
interactions of up to 6 neighbors and being a sum of roughly 100 terms (each being an 
infinite sum) . A good guess is that the explicit form of this hamiltonian will never appear 
explicitly in the literature."^ 

It is interesting to investigate the general circumstances permitting the existence of a 
ladder operator in integrable systems. For classical integrable systems, the most common 
recursive structure is the one related to their bi-hamiltonian character (see e.g. [22]). Recall 
that the existence of two distinct hamiltonian structures is an almost universal signature 
of integrability is for classical hamiltonian models. Take for instance a continuous system 
described by a set of fields collectively denoted by (p. The bi-hamiltonian character of the 
system translates into the equalities 

dap = &, H} (1) =Pij^ = {<P, H'}{2) = P 2 -^-, (1-9) 



The three-spin charge Q3 for the XYZ model was calculated in [9]. The explicit form of the 
four-spin charge appears in [19]. For the isotropic (XXX) model it can be also obtained, as a 
special case, from the four-point charge of the t-J model, written down in [20]. 

But see section (5.3)! 
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where H : H' are two distinct hamiltonians, and p (i = 1,2) denotes the differential oper- 
ator, defining the equal-time Poisson bracket { , }($) via 

Mx), <p(y)} (i) = P t (x)S(x - y). (1.10) 

If H' has higher order than H in terms of an appropriate grading, it follows that the charge 
H n+ i can be calculated from the lower order one H n by 

pr 1 ^^ = ^. (1.11) 

dip dip 

The bi-hamiltonian structure furnishes then a recursion operator, P 1 _1 P2, that acts, modulo 
functional derivatives, as a sort of ladder operator on the conserved charges. 

Constructing the quantum version of a classical integrable hamiltonian system 
amounts to a quantization of either of the two Poisson structures. The point is that 
the quantization of the system via Pi gives a theory which, in general is quite distinct 
from the one obtained from the quantization of P%. Hence, at the quantum level, the 
bi-hamiltonian character is irremediably lost.Q As a consequence, a priori there are no rea- 
sons to expect a recursive scheme for the calculation of the conserved charges in quantum 
integrable systems. 

As we just argued, the ladder property (1.5) does not reflect a sort of accidental bi- 
hamiltonian structure: the action of B is clearly not a direct quantum generalization of 
(1.11). The classical counterpart of (1.5) would be 

{-B, -ff n }(i) = H n+1: (1-12) 

for some Poisson structure { , }(j) corresponding to the commutator appearing in (1.5). 
This equation has a simple, but somewhat unfamiliar interpretation in soliton theory: it 
reflects the presence of a hamiltonian master symmetry of degree 1, which we denote by r±. 
For systems whose evolution equation takes the form 

<Pt = K(<P), (1-13) 

a mastersymmetry T\ is defined by 

[[ri,K(ip)],K(ip)] = 0, (1.14) 



This has apparently first been observed in [23] in the context of the quantum Korteweg-de 
Vries equation. 



(where the commutators should be understood as Lie derivatives). A mastersymmetry is 
hamiltonian if it can be written as 

n = P^ (1.15) 

dip 

for some H n . If Pi is a fundamental hamiltonian operator, meaning that it cannot be 
factorized into a product of powers of other hamiltonian operators, the existence of a 
hamiltonian mastersymmetry of degree 1 is the exception rather than the rule.H But in 
the Landau-Lifshitz model, which corresponds to a continuous, classical version of the 
XYZ model, such a mastersymmetry does exist [24]. This, in a sense, is the reason for the 
existence of a ladder operator for the XYZ quantum chain. 

In the XYZ context, the boost property (1.8) is actually rooted in the Yang-Baxter 
equation itself, or more precisely, in the fact that the transfer matrix is a product of the 
R matrices, satisfying the Yang-Baxter equation [17]. In other words, the Lax operator 
and the R matrix are identical. In the terminology of [8], such models are said to be 
fundamental. Somewhat paradoxically, that relation between the boost operator and the 
Yang-Baxter equation might suggest that the existence of a quantum ladder operator is 
a general feature of integrable quantum chains! However, since not all integrable spin 
models are fundamental, the existence of a ladder operator B satisfying (1.5) cannot be 
the hallmark of quantum integrable chains. In particular, no such operator exists for the 
Hubbard model (as shown below); this model is indeed not fundamental in the above sense. 
There are other examples of integrable chains without ladder operator. For instance, for the 
trigonometric (or hyperbolic) long-range interaction version of the XXX model presented 
by Haldane and Shastry [25], no such operator is knownB 

It is of great interest to study the continuous limit of the ladder operator B. In- 
deed, in appropriate sectors, ferromagnetic or antiferromagnetic, the XYZ model reduces 
either to the quantum nonlinear Schrodinger equation or the Thirring model (equivalent 
to the sine-Gordon equation by bosonization, itself related to the quantum Korteweg-de 
Vries equation by a Miura-Feigen-Fuchs transformation). Although the quantum inverse 



This statement will be made more precise in appendix A, where the concept of mastersym- 
metry is also briefly reviewed. At this point we simply mention that (1.15) is not verified for most 
bi-hamiltonian systems if Pi a fundamental hamiltonian operator. 

The original argument establishing the existence of B for the XYZ model, which is reviewed 
in section (2), breaks down for a chain with long-range interactions. 



scattering method furnishes a systematic method for calculating the conserved charges of 
these theories, it is again rather impractical. Hence, the existence of B in the XYZ model 
gives a hope that a ladder operator could exist also in these continuous models. However, 
as we will show, this hope does not materialize. 

1.2. An overview of the present results 

Our aim in this work is to unravel the general pattern of the conserved integrals. That 
amounts to first determining which terms can appear in the charge of a given order, and 
eventually fixing the values of their coefficients. Our strategy is rather straightforward: 
from the explicit expression of the first few conserved integrals, we identify the general 
pattern, which is subsequently established by means of a recursive argument. In some 
favorable circumstances, all the coefficients can be calculated explicitly. This occurs for 
two special cases of the XYZ model: the XXX model (X x = X y = X z in (1.1)) and the 
degenerate XY model (X z = 0). The conserved quantities in the latter case are well-known 
(the model being equivalent to a free fermion theory) and have been presented in [26, 27, 
28, 29]. For the infinite XXX model, a general expression for the conserved integrals was 
first reported (without proof) in [30] . Originally unaware of this work, we have rediscovered 
it in a modified form, presented in [31], with a sketch of the proof. In this work we give 
the details of the the proof. We also extend this construction to the isotropic, su(M) 
invariant generalizations of the XXX chain. In the general anisotropic case of the s = 1/2 
XYZ model, the functional form of the integrals of motion is the same as in the isotropic 
case. However, the coefficients become rather complicated polynomials in the coupling 
constants. We have been able to determine exactly only some of these coefficients. No 
significant simplification of the XYZ conservation laws has been found in the special XXZ 
case (X x = X y ). 

In the absence of a ladder operator, as is the case for the Hubbard model, the gen- 
eral form of the conserved law must first be guessed, with each term multiplied by an 
undetermined coefficient. These are then fixed by enforcing the commutativity with the 
hamiltonian. Fortunately, for the Hubbard model, a simple pattern soon emerges, in which 
the building blocks of the conserved densities are recognized as the XX conserved densities 
(which is not surprising since the Hubbard model is essentially two copies of the XX model 
interacting along their ^-component.) This observation significantly reduces the number of 
a priori undetermined free parameters. With the normalization we have chosen for the pa- 
rameter U (cf. eq. (1.2)), all coefficients turn out to be simply ±[/^, with £ a non-negative 
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integer. Furthermore, the different terms in H n can be grouped into classes, where all the 
terms of a given class have the same coefficient and can be obtained from each other by 
a relative lattice translation of its components; these classes are in correspondence with 
the partitions of the integers n — 2r > 0, where r is again a non-negative integer. The 
coefficients of the terms independent of U and those linear in U are found exactly.!!! 

Appendices are devoted to the study of some continuous theories related to the XYZ 
model. Appendix A is concerned with a general discussion of mastersymmetries in classical 
soliton theory, followed by a presentation of the results pertaining to the classical Landau- 
Lifshitz equation. Although this appendix is mainly composed with review material, it 
contains several novel observations. In particular, it determines the conditions under 
which classical mastersymmetries and a quantum boost operator can be related. This 
provides a way of detecting classical integrable models whose quantum version will contain 
a ladder operator. Next, by a direct analysis presented in the appendix B, we establish 
the evanescent nature of the ladder operator in the continuous limit of the XYZ model, 
by disproving the possible existence of ladder operators both for the quantum nonlinear 
Schrodinger and the quantum Korteweg-de Vries equations. 

Many of the calculations of the higher order charges in this work, which are often quite 
demanding computationally, were performed using a set of specialized "Mathematica" 
routines that we developed. 

1.3. Motivation 

What is the motivation for constructing explicitly the conserved charges of the quan- 
tum integrable chains? First of all, it forces us to look at these systems from a point of 
view different from the Bethe ansatz. New perspectives can obviously lead to a deeper 
understanding of these models and more generally of the nature of quantum integrable 
systems, for which, after all, relatively little is known. 

Another motivation is related to model building. It is always of interest to be able 
to mimic the gross features of a physical problem in terms of an integrable system. A 
simple way of producing new integrable physical models in the context of an integrable 



This implies the knowledge of the explicit form of all conservation laws of a hypothetical 
physical model that would be described by the Hubbard hamiltonian but with a parameter U 
satisfying U 2 = 0, a constraint that does not spoil the integrability of the model. Such a constraint 
can be realized e.g. by having U of the form U = a/3, with a and (3 constant Grassmannian 
parameters. 

10 



hierarchy, is to take as the basic hamiltonian, some linear combinations of the conserved 
charges, choosing the precise combination that produces the desirable dynamical aspects 
of the problem under consideration. The simplest illustration of this idea is the description 
of the XXZ model in a magnetic field. It is obtained by adding to the original hamiltonian, 
a term proportional to the lowest order conserved charge, namely the ^-component of the 
total spin, which is multiplied by a constant coefficient having the physical interpretation 
of an external magnetic field. Frahm has proposed recently a more sophisticated realization 
of this program, by considering a XXZ-type model with basic hamiltonian Q2 + cQ^ [32] . 

The knowledge of explicit formulae for the conserved charges could also be useful in 
finding exact solutions for small finite chains (see e.g. [33]). 

Our original motivation for this project was to see to what extent information on 
the conserved charges of the continuous quantum integrable systems could be extracted 
from those of their lattice regularization. In particular, if the explicit formulae for all the 
XXX conserved charges are known, does it mean that all the conserved charges of the 
corresponding continuous theory are also known? However, the continuous limit of the 
higher order charges is not straightforward. The problem is essentially that products of 
operators in the continuous theory have to be regularized in some way (e.g. via normal 
order) and no clear print of this regularization appears to be encoded in the lattice version. 
A full discussion of these questions is deferred to a subsequent work. 

2. The Yang-Baxter equation, commuting transfer matrices and the boost 
operator 

We first briefly review the relation between the conserved charges of a quantum spin 
chain with the commuting family transfer matrices and the Yang-Baxter equation. Since 
we do not consider in this work non-trivial boundary effects, the following discussion is 
confined to infinite chains or finite chains with periodic boundary conditions. 

Consider a quantum chain with some spin-like variables Sf, a = 1, ..., d, defined at site 
i and valued in some Lie algebra, which thereby fixes the defining commutation relation of 
the model. The dynamics is governed by a hamiltonian H containing only nearest-neighbor 
interaction terms. The canonical equation of motion 

f = [ Sj ,tf] (2.1) 
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is supposed to be equivalent to the compatibility condition for the linear system 

dt lK ' (2.2) 

Si+i(t) =Lj(A)Sj, 

which reads 

AT. 

-^ = U i+1 Li - LiUi. (2.3) 

U and L are d x d matrices depending upon a complex spectral parameter A. Finally, the 
Lax operator L is assumed to satisfy the following relation: 

R(X - ti)(Li(\) <g> Li(n)) = (Li(ji) <g> Li(X))R(X - n), (2.4) 

where R is a <i 2 x d 2 matrix with normalization -R(O) = I. The above tensor product lives 
in the space V\ <8> V2 where Vi ~ V2 — V is the <i-dimensional space on which the d x d 
quantum matrix operators L$(A) act. To indicate that -R acts on V\ <E> V2, we rewrite it as 
i^ViVb or ^12. Equation (2.4) implies as a compatibility condition the famous Yang-Baxter 
equation: 

R 23 (X - ^R 12 (X)R 23 ^) = R 12 {ji)R 2 z{X)R 12 (X - fi). (2.5) 

where V3 ~ V. It is usually written in terms of the variable R related to R by 

R = PR, (2.6) 

where P is a transposition operator, i.e. 

*,J jfe «,J = *Hk~ V 'J 

The relation (2.5) directly implies that the transfer matrix T(A) 

T(X) = Tt Vo L n (X)...L 1 (X) (2.8) 

commutes with T(p). Expanding the transfer matrix as 

lnT(A)=2zV — Q n+1 , (2.9) 

one obtains a family of mutually commuting quantities {Q n }'- 

[Qn,Qm}=0. (2.10) 
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The matrix elements of the operator Li (A) are functions of the quantum variables S*, 
themselves d x d matrices. Matrix entries of Li act then on a <i-dimensional space, called 
the quantum space and denoted by Vo- Because Vo — V, we can set 

L i (\) = R i0 (2.11) 

(up to a normalization factor). In that case, (2.4) becomes 

R 12 (X - ^R 10 (X)R 20 (^ = R2o(n)Rio(\)Ri2(\ - n), (2.12) 

which, thanks to (2.6), is the same as (2.5). 

Spin models for which the Lax operator is the same as the R matrix are called fun- 
damental [8]. In such cases, the quantum hamiltonian can be expressed directly in terms 
of the R matrix. Indeed, the transfer matrix is then a product of R matrices, one for each 
site: 

T(A) = Tv Vo R N0 (X)...R 10 (X), (2.13) 

so that 

dT{x) \ -Vtt p p dRl ° {x) \ IP P )P P 

1 v (2.14) 

EdRa-\-\\X) . 
^ \ x=0 Tr Vo P Nfi ...Pi,o 

% 

(in the first line we used the fact that Pf = I). Up to a constant, this is simply the 
quantum chain hamiltonian 

fl = -T(0)-'^U = -E^%^U A.m = EW (2.15) 

We next briefly review the relation between the Yang-Baxter equation and the boost 
operator [17]. The starting point is (2.12) with 1 — ► k, 2 — > A; + 1 and X = u + fx. The 
derivative of this expression with respect to v evaluated at v = yields 

; Rok+i(/J>) — Rok ; = [hk,k+i, RokRok+i > (2-16) 

d[i d\x 

where hk,k+i is defined in (2.15). This result is then multiplied from the left by 
Iln<fc R 0n{v) and from the right by Iln>fc+i Ronif*)'- 

(URoM) d -^(l[RoM)-( n RoM)^f^( n R oM) 

n<k " n>k n<k + l " n>A; + l (2-17) 

= [^ fc>fc+1 ,T(/i)]. 
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Multiplication by k followed by a summation over k from — oo to oo leads then to (1.8), 
with B defined by (1.7). From (2.9), this directly implies 

[B,Q n )=Q n+1 . (2.18) 

For finite chains with periodic boundary conditions, the argument is the same except that 
the coefficients are understood to be defined modulo N, the number of sites. 

The logarithm in (2.9) ensures the locality of the charges Q n , i.e. that Q n does not 
contain interactions between spins at distances greater then n lattice units. As shown by 
Liischer [9] for the XYZ model, these charges can be put in the form: 

Qn= Yl G£-l(*l, ■■■,*»-!), (2-19) 

{il,...,»n-l} 

where the summation is over ordered subsets {ii, . . . ,i n -i} of the chain, and G T is a 
translation covariant and totally symmetric function, with the locality property: 

G^(n,...,z n ) = 0, for \i n -i x \>n. (2.20) 

For the infinite XXX chain, further properties of the conserved charges, including their 
completeness, have been proved in [34] B 

In the rest of the paper we will reserve the notation {Q n } for the charges defined by 
(2.9). Note that a rescaling of the spectral parameter A — > a -1 A results in a multiplicative 
redefinition of the Q n 's : Qn —* & n ~ 1 Qn- For the XYZ model it is convenient to use the 
freedom to rescale the spectral parameter so that the logarithmic derivative of the transfer 
matrix becomes: 

-T- 1 (0)r(A)| A=0 = \f>2 = ^E^K+i + VK+i + A,oJoJ +1 ]. (2.21) 

j'eA 

With this normalization, the boost operator for the XYZ model reads: 

B = ^E^^i+i + A ^K+i + A ^i^i+i]- (2-22) 

j'eA 

The structure of the conservation laws is however more transparent in another basis, which 
may be obtained by taking an appropriate linear combinations of Q n 's. Such a basis will 
be denoted {H n } and chosen so that H n = ,^ 2 y_ Qn+ a linear combination of lower order 

Qi's. 



To our knowledge, no such completeness proof exists for the anisotropic case. 
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3. The conserved charges of the XYZ model: notation and generalities 

3.1. Notation 

We begin by introducing the necessary notation. We consider either a finite spin 
lattice with periodic boundary conditions (A = {1, . . . , iV}, with N + 1 = 1) or an infinite 
one ( A = Z) . Our constructions applies equally well to both cases, provided that addition 
in A is understood modulo N for the finite chain. We now introduce two general objects 
that enter naturally in the construction of the conserved charges: clusters and patterns. 

A sequence of n lattice sites C = {i±, ...,i n }, with i\ < i 2 < ... < i n , will be called a 
cluster of order n. A disordered cluster will refer to a sequence with the ordering condition 
relaxedB A cluster starting at i\ and ending at i n , and containing non-adjacent sites has 
k = i n — ii + 1 — n holes (sites in between i\ to i n that are not included in C); k = for a 
cluster containing only contiguous spins. The set of all clusters of A of order n with k holes 
will be denoted as C^ k \ For instance, C^V contains {1,2,3,6}, {1,2,5,6}, {1,4,5,6}, 
{1, 2, 4, 6}, {1, 3, 4, 6}, {1, 3, 5, 6} and all their translations. 

Let S( n > denote the set of all sequences of n basis spin matrices of the form 

£ = {<,...,<;}, (3.1) 

where Oj = x, y, z, and i\ < i 2 < ... < i n - Clearly, such a sequence is completely specified 
by the cluster {ii, . . . , i n } and by the assignment of group indices, {ai, . . . , a n }, which 
we call a pattern. We denote by V^ n > the set of all n-spin patterns, i.e. sequences of n 
elements taking values in the set {x, y, z}. There exist natural projections, assigning to a 
spin-sequence S G S^ n ' its corresponding pattern or cluster: 

p : S in) -»■ V {n) and p' : S {n) -> \JC {n ' k) (3.2) 

k 

given by: 

p(S) = {ai,a 2 ,...,a n }, p'(S) = {ii,i 2 , ■ ■ ■ ,i n }- (3-3) 

We will write C 71 " to denote the sequence S with cluster C and pattern tv. 



This definition of a cluster differs slightly from the one that we used previously in [31]. 
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Next we construct the n-linear spin polynomials f n , denned on basis sequences in 
<S( n ) 5 which will be the building blocks of our construction. First we define two different 
rescaling of the spin variables:!^ 



(3.4) 



^-i / \ _a 






In these variables, the XYZ hamiltonian reads: 

jeA 

Then, for any sequence of spins S G S^ n \ we define V m («S) = {V^, V^, V^}(«S) (for m < n) 
as the vector product of the first m spins of the sequence, where the leftmost spin factor 
appears with a hat and all others with a tilde, with parentheses nested toward the left, 

Vi(S) = & h , 

V 2 (S) = a ix x a i2 , 

V 3 (5) = (& it x a- i2 ) x a i3 , (3.6) 

V m (S) = V m _i(<S) x a im . 
We also denote by a m (iT) the " direction" of V m (C 7r ), i.e: 

am(7r) = (• • • (oi x a 2 ) x a 3 ) . . . x a m ). (3.7) 

Then, we construct scalar n-linear polynomials, defined in the space of spin sequences, 
from the scalar product of V n _i(5) and &i n , multiplied by an appropriate constant g(S): 

f n (S)=g(S)(V n _ 1 (S)-a in ), (3.8) 

with fo(<S) = fi(<S) = 0. The coefficient g(S) is a polynomial in A x , X yi \ Z1 determined by 
the pattern and the positions of the holes in the cluster: 



Although these rescalings involve square roots of the coupling constants, the resulting ex- 
pressions for charges will contain only integer powers of the couplings. 
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(i)- If S contains only contiguous spins, g(S) = 1. In this case 

f n (S) = ((... ((*£ x a-) x?J)x.. OC-t) ^ <" 



j- a l ^ a 2 ,^- a 3 Z- a n-l p_a 7 






~ "n-1 ^ 

(3.9) 



where 

e(7r) = sgn[a n _i(7r)]5 a „_ lU))a „ (3.10) 

with the convention sgn(0) = 0. Hence, e(7r) can be either or ±1. For example, we have: 



(3.11) 



/2({o-f,crf +1 }) = A x ofof +1 , 

/ 3 ({<,crf +1 ,< +2 }) = A x A y ofof +1 of+2, 
/ 3 (K,af +1 ,ar +2 }) = 0, 

/ 4 ({<,af +1 ,af +2 ,< +3 }) = -A X 2 A, CTfcTf +1 (7f +2 (7f + 3. 

(ii)- For clusters containing holes, the coefficient g(S) is given by 

g(S)= n A * A * A ' , (3.12) 

where L(j) is the number of sites in C to the left of site j and cll(j){^) has been defined 
above in (3.7). For example, 

; 3 (K,af +2 ,exf +3 }) = (A,A,)A x A y ofafaafa, (3.13) 

/ 3 (K,af +1 ,af +3 ,af +4 }) = -(A.A^A^A, afaja^af, 

(where the coefficient g(S) is indicated in parentheses). 

Finally, we define the linear spaces T n ,k spanned by the translation invariant local 
polynomials built from / n 's:&il 

Fn, k = Yl «( ff ) E /"C^)' ( 3 - 14 ) 

7rE p(n) cec™> fc 

where a(7r) are some real-valued coefficients. 



J- n ,k is a linear subspace of dimension 3 • 2 n 2 of the space of observables; since the product 
of two local polynomials is not local, it is not an algebra. 
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In the isotropic case, i.e. for the XXX model, all the coupling dependent coefficients 
become one, which results in a significant simplification. We can then define the polyno- 
mials f n whose arguments are clusters in C n ' fc , as follows: 

fn(C)= J2 ~fn(<^)- ( 3 - 15 ) 

7reP(") 

In particular, omitting the argument of the / n 's for brevity, we have 

fo = fi = 0, 

H = Ci 1 -Ci 2 i 

/3 = (o*ii x <7i 2 )- a 'i3i 

U = (foi x 0"i 2 ) X O'^V 



(3.16) 



The multilinear polynomials / n 's and / n 's satisfy an interesting property which is that 
the dot product can be placed at an arbitrary position, provided that parentheses to its 
left (right) are nested toward the left (right), e.g: 



h = (((o"ii x (T i2 ) x a i3 ) x o- i4 )-er is = ({a- il x cr i2 ) x a i3 ) ■ {a u x <r ib \ 
= a ix ■ (er i2 x (<7; 3 x (cr i4 x <T is ))). 



(3.17) 



This is a direct consequence of the familiar vector identity: 

(AxB)-C = A-(B x C). (3.18) 

3.2. The structure of the XYZ conservation laws 

The conserved charges of the XYZ model have the general form: 

[n/2]-l k 
Qn= Yl E^-2 M , (3-19) 

k=o e=o 

where the square bracket stands for integer part, and F™_ 2k e is an element of T n -ik,(.- 
More precisely, 

K- 2k , e = J2 e n-2k,M E fn-ik(n, (3.20) 
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where 0™_ 2k e is a polynomial of degree 2(k — I) in the coupling constants A x , X y , X z . It is 
convenient to think of Q n in terms of a triangle, consisting of the operators F™_ 2k e : 











77171 








^n-2,1 


pn 
^71-2,0 






r n-4,2 


^n-4,1 


r n-4,0 




r n-6,3 


fpn 
r n-6,2 


r n-6,l 


pn 
^71-6,0 


T-in 
r n-8, 4 


r n-8, 3 


fpn 
r n-8,2 


-^n-8,1 


77171 
^71-8,0 



(3.21) 



The bottom edge of the triangle is formed by the sequence {.F^Jk-o,.. r n /2]-i if n is even, 
or {i 71 ^j t }/ i;= o,...,[n/2]-i if n is odd. Note that linear combinations of the form 

[n/2]-l 
-H" n = 2^ Pn-2fc(A x , A y , \ z )Qn-2k, (3.22) 

fc=0 

where p n -2k is an arbitrary polynomial of degree n — 2k in the coupling constants preserve 
the form of the triangle, redefining only the polynomials 0{tx). 

The problem is how to determine the coefficients 9{ix). In section (5), where we 
prove ( |3.19| ), we give a recursive method to calculate F™_ 2k t . Some of the coefficients 
are easy to calculate: e.g. one finds that 9™_ 2k k (jc) = (n — 2)! for all patterns tt (thus 
F n-2k.k ~ Z)7re-p(™- 2fe ) Z)c6C(™- 2fe . fe ) fni^))- This fixes the coefficients on the left edge of 
the triangle. Unfortunately, finding all the coefficients in the general anisotropic XYZ case 
is a very difficult task and we have not succeeded in disentangling the equations for the #'s 
apart for those on the left edge. However, there are two particular cases in which all the 
coefficients are relatively simple and can be determined exactly, namely when X z = or in 
the isotropic limit where A x = X y = A 2 . In the XY model, which is analyzed in section (6), 
the triangle degenerates into its right edge. On the other hand, for the XXX model, treated 
in detail in the next section, all the coupling-dependent coefficients become constants. In 
that case, by an appropriate transformation of the type Q3.22 ), one can eliminate all the 



terms on the right edge, except for F™ . This decoupling provides enough simplification 
in the recursion relations for the coefficients to allow for an exact solution. 

For the XXZ model, the recursion formulae are somewhat simplified, being dependent 
only on a single parameter X z /X x . However, this simplification is not as dramatic as in the 
XXX case, and it does not seem to permit a solution in closed form for all the coefficients. 
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4. The explicit form of the conservation laws in the XXX case 

4-1. The general formula 

For the isotropic XXX model the only quantities that enter in the construction of the 
local integrals of motion are 

F n , k = J2 ^(C). (4.1) 

Cec(™> fe ) 

Quite remarkably, the linear combination of -F n ,fc's that builds up H n is encoded in a simple 
tree, in which the vertices are labeled by the F n ^'s. This tree is displayed in Fig. 1. 



Fn-4,2 



Fn-4,l 




• • • • 



Fig. 1. The tree structure corresponding to H n . The tree stops with the terms F2/ 
(F3/) when n is even (odd). Note that the trees corresponding to H 2 m and H2 m +i have 
identical structure, differing only in labeling of their vertices. 



The sum of all the vertices of the tree, with each vertex contributing with unit weight, 
gives H n . Summing up all the vertices of the H n tree leads then to the following expression: 



[n/2]-l k 

H n = F n fi + 2_^ / v <*k,lF n 

k=l £=1 
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-2kl, 



(4.2) 



where the coefficients ot.k,t are defined recursively by the relation: 

k 

m=e-i 

with ckii = 1 and o^o = 0. Notice that au,x = Oik,2 f° r k > 2. These coefficients turn 
out to be generalized Catalan numbers, as can be seen by rewriting the recurrence relation 
(fO|) in the form: 

otk,t = Oik-i,l-i + c*k,i+i, (4.4) 

with the convention a.k,i = if £ > k. This is the defining relation for the generalized 
Catalan numbers, cek,£ = C2k-i-i,i> with C n , m given by 

where I , J are the binomial coefficients, with p = \{n — m + l)/2], m + n odd and 

m < n+2. In particular, ak,i = C*2fe,i reproduce the familiar sequence of Catalan numbers: 
{1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862 . . .}. The tree in Fig. 1 is known as a Catalan tree. 
The first few conserved charges are: 

Hi = ^2,0, 
H3 = i*3,0j 

H 4 = F 4)0 + F 2jl , 
H5 = F 5)0 + -F 3; i, 

Hq = F 6) o + -^4,1 + ^2,2 + -^2,1 5 

H7 = -F 7)0 + -F 5jl + -^3,2 + ^3,1) 

Hs = -^8,0 + -?6,1 + F&,2 + -?4,1 + -^2,3 + 2^2,2 + 2^2,1, 

Hq = -£9,0 + ^7,1 + -^5,1 + -^5,2 + -^3,3 + 2-F 3i 2 + 2F 3]1 , 

Hiq = -Fio.O + -^8,1 + -^6,2 +-^6,1 + -^4,3 + 2-^4,2 + 2F 4jl 

+ F 2 ,4 + 3F 2 ,3 + 5F2,2 + 5F 2 ,i. 



(4.6) 



We note that O can be put in the form ( gT9|) with F™ = F n , , F£_ 2kfi = ° for 
< k < [n/2] - 1, i^_ 2M = a M F n _ 2M . 

The above construction gives iV — 1 charges {-ff2 5 • • • -Hjv } f° r the XXX chain of length 
iV with periodic boundary conditions. We can add to the set {H n } any of the three 
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components of the total spin, Hf = J2jeA^j- As the ^n,fc' s are invariant under global 
spin rotation, [iff, iJ n ] = 0.t3 For the infinite isotropic Heisenberg chain we have then an 
infinite sequence of independent local charges {H® , H 2 , ..., H n , ...}. 

4-2. Details of the proof 

We now present the direct algebraic proof of the mutual commutativity of all the 
{H n Ys. We will first show that [H 2, H n ] = 0, by evaluating the commutator: 

[H 2 J n (C)}=J2[<rr<rj+iJn(C)}. (4.7) 

jeA 

The proof is rather lengthy and some of its steps are quite tedious. It is divided into nine 
parts. 

a) Three useful identities 

To proceed, we need a few simple identities for commutators in the su(2) tensor 
product space. Let A = {a^ , a y -, a z A and B = {a^, o v k , a^} denote the basis of two different 
(j 7^ k) su(2) factors, and let L = {L 1 , L 2 , L 3 } and R = {R 1 , -R 2 , R 3 } be arbitrary tensor 
products commuting with all components of both A and B. The following identities hold: 

[A • B, (B x L) • R] = -2z((A x B) x L) • R, (4.8) 

[A ■ B, (A x B) • R] = 4z(B • R) - 4z(A • R), (4.9) 

[A • B, ((L x A) x B) ■ R] = 2z(L x B) ■ R - 2z(L x A) • R. (4.10) 

These identities are simple consequences of the basic su{2) multiplication rule for the 
components of A, B: 

A a A b = S ab + ie abc A c . (4.11) 
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Thus the charges (^2) survive in the presence of a magnetic field coupled to any of the 



components of the total spin. 
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b) Commuting (Xj<Xj+i with f n (C) 

The above identities are helpful in evaluating the commutators of the basic links 
(TiCTi+i with f n (C). It is remarkable that all such commutators are of the form f n ±i(C), 
where the cluster C can be obtained by applying a number of simple transformations to 
the cluster C = {ii, . . . , i n }: 

ii-lC = {h ~ l,^i, - • -,i n }, 

(4 - 12) 

C-r- = {ii, . . . , ik-i, ik+i, ■ ■ ■ , in}, 

^ij — >ifc ii = \*1 > ■ • • ) *J — 1 5 \*/c» *^y ) "J+lj ■ ■ • ) <"nf- 

Note that the last operation is defined only if i/-, it are not in C^; the brackets indicate the 
change in nesting of vector products when the expression appears as an argument of the 
polynomial f n +i- As a result, this fourth operation produces a disordered cluster, unless 
ij is at the extremity, i.e. j — 1 or j — n. Although this may not be completely clear at 
this point, it will be illustrated shortly by an example. On the other hand, the first three 
transformations always yield ordered clusters. 

For C E c( n ' k > the terms arising in (|4.7|) correspond to clusters with (??/, k') — (n + 
1, /c), (n + 1, k — 1), (n — 1, /c), (n — 1, k — 1). More precisely, using the above identities, 
the different terms in ( f4.7|) are found to be: 



(4.13) 



[o*i-l • (Ti x , fn(C)] = —2ifn+l(ii-lC), 

[cri n -i ■ cr in ,f n (C)] = 2if n+ i(C in+ i), 

[0*1 ■ 0*1+1, fn(C)] = 4i/ n _i(C-) - 4z/ n _i(Cr^ ri ) if i 2 = i\ + 1 

= 2if n+1 ( ilil+1 C~) if J 2 ^i + 1, 
[<Ti n -i-<Ti n ,fn(C)]=4if n - 1 (C r -~~)-4if n - 1 (C^) if i n -i = i n - 1 

= -2z/ n+ i(C-i„_ H „) if i„_i ^i„ - 1, 
[<Ti 3 ■ 0*,-+i, f n (C)] = 2if n - 1 (C^) - 2if n - 1 (C- — ) if ij,ij+i G C 

= -2if n j r i(C ij ^i j i j j r i) if ij eC, ,ij+i£C, 

where in the last case, it is assumed that j ^ 1, n — 1. 

The last term in 4.13 contains contributions from disordered clusters. For example, 
consider the cluster C = {1, 2, 4, 5} and the commutator of 

f 4 (C) = cr 1 ■ (cr 2 x (cr 4 x cr 5 )) = (er 2 x (cr 4 x <x 5 )) • <n. (4.14) 
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with er 2 • cr 3 . Taking for f± the second expression above allows us to use (|4.8|) directly, 
with the result 

[cr 2 ■ cr 3 , (<x 2 x (cr 4 x cr 5 )) ■ o-j] = -2z((cr 3 x cr 2 ) x (cr 4 x cr 5 )) • cri (4.15) 

We see clearly that the nesting of the spins has been modified (that is the parentheses 
cannot be all nested toward one side or the other). In other words, the result cannot be 
written in terms of /s({l, 2, 3,4, 5}). As this calculation illustrates, the special property 
( p,18|) makes the direct application of the identities (|4.8| ) -( |4.10|) straightforward. 

c) Summation over links inside a cluster 
Adding up all these terms gives: 

[if 2 , fn(C)] = a n+1 , k (C) + b n -i,k+i(C) + d n+ i jfc _i(C) + e„_i,fe(C) + r{C), (4.16) 

where the various quantities appearing in this expression are defined by: 

a n +i,k(C) = —2if n+ i( il -iC) +2if n+ i(Ci n +\), 
6 n _i jfc+ i(C) = -4i/ n _i(C£)<S il +i, ia +4i/ n _i(Cr— ;)^ n _ 1 +i,i„ 



+ 2i X)[/»-i(^) - /»-i(Cj^)]Vi,*i+i' 

J =2 
^n+l.fe-l(C) = 2z/ n _|_i(j 1 j 1 _|_iC^)(l — (^ j 2 _i) 

- 2z/ n+1 (C r ,„_H„)(l - 6 in , in _ 1+1 ), (4.17) 

e n -i,fc(C) =4zf n _i(C-)5 il+ i j i 2 -4z/ n _i(C-)5 i „_ 1+ i )i „, 

n-2 

r(C) = -% Y^ fu+i{.C ij ^ jij+1 )(l - 6i. +lti . +1 ). 

J =2 

n-1 
+ 2z 2_^ /n+l (Cij ^ij -lij ) (1 - ^_i,ij-l) 

i=3 
Caution: in the above formulae, one should distinguish carefully between sites ij + 1 and 
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d) Summation over clusters 

With the exception of r(C), all of the terms on the right hand-side of (|4.16|) involve 
only regular nested product of the full set of spin variables. Fortunately, when summing 
up over all clusters in C^ n,k \ the unwanted contributions from improperly-nested spin 
sequences cancel: 

E KC) = 0. (4.18) 

This can be seen as follows, r-type terms arise only from clusters containing holes. Consider 
a particular cluster C having a hole at site j but not at site j + 1. Let us write its associated 
polynomial / as 

/(C') = (Lx,t j+1 )-R, (4.19) 

where L stands for the contribution of all the spins at the left of the site j and R for that 
of the spins at the right of j + 1. The contribution of this cluster to the commutator with 

(TjCTj^i IS 

[o-j(Tj + i, (L x cr j+1 ) ■ R] = — 2z(L x ((Tj x (Tj+i)) ■ R. (4.20) 

In the summation over all possible clusters having the same number of spins and same 
number of holes, we will encounter another cluster C" , which differs from C only in that 
the hole which was in position j now appears in position j + 1. The / polynomial of this 
cluster being 

f(C") = (L x aj ) ■ R, (4.21) 

its contribution to the commutator with (Tj(Tj + i is 

[(Tj<Tj +1 , (L x o-j) ■ R] = -2z(L x (cr j+1 x (Tj)) ■ R, (4.22) 

which is exactly the negative of (4.20). Hence, in the summation over all clusters of a fixed 
type, all the contributing r-terms cancel two by two. 

It then follows that the commutator [H2,F n ^] can be expressed as a linear combi- 
nation of terms of the type F n i ,*/ . Equivalently, we may interpret this commutator as a 
transformation of the set of all clusters into itself. Now, with an argument similar to the 
one used to show the cancellation of the r terms, we can prove that the sum over the e-type 
terms also adds up to zero: 

J2 e n _ ljfc (C)=0. (4.23) 

Cec( n ' k i 
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Let us illustrate this calculation by a detailed example. Take iV = 7, n = 4 and k = 1. 
The different clusters C^' 1 ' are 



{C^ l) } N=7 ={{1, 2, 4, 5}, {2, 3, 5, 6}, {3, 4, 6, 7}, {1, 3, 4, 5}, {2, 4, 5, 6}, 
{3, 5, 6, 7}, {1, 2, 3, 5}, {2, 3, 4, 6}, {3, 4, 5, 7}}. 

The e 3) i terms for the first three clusters are 

4i[/ 3 ({2, 4, 5})-/ 3 ({l, 2, 4}) + / 3 ({3, 5, 6}) - / 3 ({2, 3, 5}) 
+ / 3 ({4,6,7})-/ 3 ({3,4,6})], 

those for the following three clusters are 



where 



A n +l,k — 2_^/ a n+l,k(C), 

C e C(™,fe) 

B n -l : k+l = 2_^t ^n-l,k+l(C)i 

C e C(™,fe) 

D n +i,k-l = 2_^ dn+i,k+i(C)- 

C e C(™,fe) 



(4.24) 



(4.25) 



4i[-/ 3 ({l, 3, 4}) - / 3 ({2, 4, 5}) - / 3 ({3, 5, 6})], (4.26) 

and for the last three, these are 

4i[/ 3 ({2, 3, 5}) + / 3 ({3, 4, 6}) + / 3 ({4, 5, 7})]. (4.27) 

Adding up all these contributions, we see that terms cancel in pairs. 
Therefore, 

[if 2 , F n> k] = A n+ i : k + -Dn+i,fc-i + B n -i t k+i, (4.28) 



(4.29) 



e) The special case where n = N 

In the case where n = N, the calculation of ( |4.7|) is slightly different, as a new type 



of term appears in the commutator of the hamiltonian and the highest order term in Hn, 
i.e. 

iV 

F N ,o = ^2f N {Cj), (4.30) 

i=i 

where 

C j = {j,j + l...J + N-l}. (4.31) 
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This commutator reads: 

[H 2 , f N (Cj)] = 9j + e N - hl (Cj) + b N - ltl {Cj), (4.32) 

where 

9j = [o-j+N-i ■ (TjJ N (Cj)] 

4.33 

— n^ka b Q b a 1 b 2 b N - 3 a N ^ 2 k f_a ai n- aN - 2 — n- a ° rr ai n- aN - 2 "l 

— Zlt t ■ ■ ■ t \ a j + l a j+2 ■ ■ ■ a j+N-l a j a j+l • ■ ■ a j+N-2)- 

It is obvious that 

JV 

Y, 9j =0. (4.34) 

3=1 



Summing up over all Cj in ( [4.32 ) gives then 

[H 2 ,F Nfi ]=B N . 1A . (4.35) 

f) Conditions for the commutativity of H n with Hi 

A simple consequence of (|4.28|) and ( [4.35| ), is that the commutation with the hamil- 



tonian changes the order of the clusters in F n & from n to n ± 1. Therefore, to prove that 



the charge H n of the form (|4.2| ) commutes with the hamiltonian, it is sufficient to show, 
for any n < N, that: 

(i) [i?2, -fn,o] does not contain terms of order n + 1, 
(ii) [H 2 , F n ,k + ^2i=i F n -2 t i] does not contain terms of order n — 1. 

The first assertion is easily verified for n < N since A n _|_i i o = 0, by translational symmetry, 
that is 

A n+lfi = 2iJ2i-fn+i({J - 1, J, ...,j + n}) + f n+ i({jj + 1, -, J + n + 1})] = 0. (4.36) 
jeA 

For n = N, (i) immediately follows from ( [4.35| ). The second assertion is more difficult to 
prove. It is equivalent to: 

fc+i 
S n _i >fc+1 + J2( A n-i,e + Ai-M-i) = 0. (4.37) 

As the sum above contains terms corresponding to clusters of order n—1, with hole numbers 
ranging from to k + 1 , it can vanish if and only if all terms with a given number of holes 
cancel, i.e. 

B n -i t k+i + Ai-i,fc+i = 0, (4.38) 
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A n - lt i + D n _ ltt = 0, (l<£<k) (4.39) 

fln-1,0 = 0. (4.40) 



Q4.39|) and ( |4.40|) are direct consequences of (4.17). In fact, they follow immediately once 



we notice that D n -i^ can be rewritten as 

D n _ M = - J2 a n+iAC) (4-41) 

cec(™> fe ) 

The proof of ( |4.38| ) is much more tedious. 

g) A resummation trick 

In order to demonstrate ( f4.38|) , one first has to rewrite the B term in a way that can 
be compared with the form of A. B naturally decomposes into two parts; 

B n _ hk+l = B n _ lk+1 + B n _ lk+1 , (4.42) 



with 



B ( n \ k+1 = -4i J2 lfn-i(C r2 )S h+ i, i2 - /n-i^M^+MJ, (4.43) 

C e CO-i,fc+i) 



and 



n-2 

,(2) 



B { n\ k+ i = 2z £ Ef/n-i(Cr) - /-iiywi- ( 4 - 44 ) 



*3 + l 



Both parts vanish identically for n = 2, 3. For n > 4, the set of clusters C^ n 1 ' fc + 1 ) can be 
decomposed into the following four classes: 



Ci = {ii,ii+i,Ko,i n -l,i n } ~ ••K * 

C2 = {h,ii+i,K ,i n _i<i n -l,i n } ~ ••if 1 o 

C 3 = {h,i2>ii + l,Ko,i n -l,in} ~ • o if{ • 

C 4 = {ii,t2 >h + l, K Q , i n -x <i n -l,i n } ~ »oif 2 o 



(4.45) 



Here -Kp (with or without prime) stands for an element of C^ n ~ 5+P ' fc + 1_ p). On the rhs we 
have introduced a symbolic description of the cluster in which the first two entries refer to 
the sites i\ and i% + 1 (• means that the corresponding site is included in the cluster and o 
means that it is not) and the last two give the information relative to the sites i n — 1 and 
i n . Everything in between is specified by the Ki's. Note that the clusters of the type C4 
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do not contribute to B n _ x k+1 due to the constraint encoded in the delta function. With 
the above schematic description for the clusters, we have: 

K p eC<"- 5 +P'*+ 1 -!>) (4.46) 

+/n-l(« o Ki o •) - / n _i(« o K( o •)]. 

Since the summation of terms containing K\ and -ftT{ runs over identical sets, their contri- 
butions cancel out, leaving us with 

B { n \ k+1 = -4i Yl [/»-!(• ° *o ••)-/»-!(• «tfoo.)]. (4.47) 

(2) 

We now turn to B n _ x k+1 . The crucial step in the present argument is the rewriting of 
(4.44) in the form 

n-2 n-2 

B^, k+1 = 2i[Pr Yl Y.fn-i{C t )-Pr £ £/n-i(C— )]• (4.48) 

CgC("- 1 ' fc + 1 ) J=2 C6C(™- 1 ' fe + 1 ) i=2 

Notice first that the constraint ij+i = ij + 1 has been relaxed. The symbol Pr indicates 
a projection operation whose action amounts to rescale the multiplicity of a given cluster 
to one. The rationale for this projection is as follows. In (4.44), the clusters C^ for a 
fixed j are all distinct due to the constraint ij+i = ij + 1. But when the constraint is 
relaxed, this is no longer so and it is then necessary to rule out repetitions. Here is an 
example. Translational invariance being not an issue here, consider a subset of C^ n,k \ 
consisting of clusters with fixed extremities (and this is what the tilde denotes below). Let 
n = 4, k = 1, i\ = 1, Z4 = 5: 

{C} ={{1,2, 3, 5}, {1,2,4, 5}, {1,3, 4, 5}}, 
{<%}<,=«,+! ={{1, 3, 5}, {1, 4, 5}}, (4.49) 

{^} 13=J2+1 ={{1,2,5},{1,3,5}}, 



1 3 ' 



so that 



£[/3(<%) - /3(<%)]<Wi = / 3 ({1, 4, 5}) - / 3 ({1, 2, 5}). (4.50) 

c 
On the other hand, since 



{C r J ={{1,3,5}, {1,4,5}, {1,4,5}}, 



?2 



{C r J ={{1,2, 5}, {1,2, 5}, {1,3, 5}}, 
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(4.51) 



we see that 

£>(<%) - MCq)] = 2[/ 3 ({l, 4, 5}) - / 3 ({1, 2, 5})]. (4.52) 

c 
By working out another simple example, with n = 5, k = 2, the reader will convince 
himself that the projection has to be done for each sum separately before evaluating their 
difference. 

We now prove the equivalence between (4.44) and (4.48). Clearly (4.48) incorporates 
all the terms in (4.44), which are those associated to clusters containing both ij and ij ■, + 1. 
The extra terms in (4.48) turn out to cancel two by two: the contribution of a positive term, 
with ij G C and ij + 1 ^ C is canceled by the contribution of a negative term associated to 
a cluster C which differs from C only in that ij + 1 G C and ij G" C . Schematically, if 

C ~ K • oK' (• at site j), 
C ~ K o •K' (• at site j + 1), 
then 



(4.53) 



C- = Ci—. (4.54) 

Having established (4.48), we now rewrite it in the form 

*iVH! = 2i[Pr £ fn-iOfy-Pr £ /»-i(fe )]• ( 4 - 55 ) 

cgc(n-i,fc+i) cec( ii ^ 1 - fe + 1 ) 

The sum over clusters can again be separated into sums over the four classes (4.45), with 
the result: 

B%l hk+1 =2iPr Yl [/»-!(• ° K • •) + /«-i(« ^i°») 



X ie C(«-5+i,fe + l-i) 



+/„_!(• o K[ f . .) + /„_!(• o K 2 ~ o 

6 2 fc 2 



-2zPr £ [/ n _ 1 (..Koo.) + / n _ 1 (..K 1 , 

x i eC(™- 5 + i ' fe + 1 - i ) 

+/n-l(» °^l°«) + /n-l(» ° -#2 



(4.56) 



ln-1 

The second and the seventh terms cancel and similarly for the fourth and the eight ones. 
The contribution of the third term being already included in that of the first one, it can be 
ignored (recall that due to the projection factor, multiplicities are irrelevant). Similarly, 

(2) 

the sixth term can be dropped. With these simplifications, B n _ x k+1 take the form 

B { *l ljk+1 = 2i J2 [/n-i(«oK ..)-/ n _ 1 (..K o.)], (4.57) 

K eC("- 5 ' k + 1 ) 
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the projection operator being now unnecessary. Adding (4.47) and (4.57) yields 

B n _ hk+1 = -2i J2 [/n-i(«oK ..)-/ n _ 1 (..K o.)]. (4.58) 

K eC(™- 5 . fe + 1 ) 

The trivial modification 

B n _ 1:k+1 = -2i ^ [fn-i(»°K »») + / n _i(««K_i ••) 

Ko,K-i (4.59) 

-/n-l(« • -f^O ° •) - /n-l(» • ^-1 • •)], 

allows us to reach the form 

B n -i tk +l = - 2z ^ [fn-l(C •) - /n-l(» C)] 

C e C(«-i,fe+i) 

= -2i £ [A-i^Wn-iC^C)] (4.60) 

C £ C(«-i,fc+i) 
= — ^n-l,fe+l) 

which thereby proves the desired result. 

h) Mutual commutativity of the H n 's and relation with the Q n basis 

Having established the commutativity of the H n J s with the hamiltonian, the final step 
in the proof of integrability is to show that all these conservation laws commute among 
themselves. For this, we will proceed indirectly by first expressing the charges H n in the 
{Qn} basis defined recursively by 

[B,Q n ] = Q n+1 (4.61) 

with Q2 — H21 and then use the fact that the Q n mutually commute, being directly related 
to the commuting transfer matrices. We first evaluate the commutator of H n (n > 2) with 
the boost operator. The result is: 

max(l,[n/2]-l) 

[B,H n ]= J2 /#°#»+i-2fc, (4.62) 

fc=0 

where. 

/#° =n - 1, 

/3 (n>2) =5 _ ^ ^g^ 

/ 3 i<£<[n/2] =~(n-2£- l)a e ,i 
31 



The simplest way of obtaining this formula is to compare the coefficients of the terms 
F n _2/ c +i j o on both sides. With this result, it is simple to obtain the expression for the 
charges Q n in terms of the {if n }'s. It is clear from fl4.62p that for even (odd) n, Q n can 
be expressed as a linear combination of the H m with even (odd) m < n: 

[n/2]-l 

Qn= J2 -yi> n) H n - 2p . (4.64) 

The coefficients 7 can be determined from a simple recurrence relation:!!^ 

1^= E 7^^" 2p) , (4.65) 

p-\-T7l = £ 

(2) 

with 7p = 5 p fi. For example, up to additive constants, 

Q 4 = 2H4, - 4ff 2 , 

Q 5 = 6# 5 - 18^3, (4.66) 

Q G = 24H G - 96H 4 + 72H 2 . 

Therefore, we see that {H n } represents a basis of the space of the conservation laws of 
the XXX chain, which can be obtained taking linear combinations of {<3n}'s. Mutual 
commutativity of H^s now follows trivially from that of the Q n 's. 

i) An alternative finale 

It is also possible to prove the involutive nature of H n without using the transfer 

matrix formalism at all, by applying an inductive argument using ( |4.62| ) and the Jacobi 

identity. We sketch this argument below. First, we note that ( |4.62[ ) can be rewritten in 

the form 

H n+1 = —!— [B, H n ] + Rn, (4.67) 

(n-1) 

where R n is a linear combination of the charges H m<n . Then, assuming that [H n , H m ] = 

for all n, m < no, we prove that [H no +i, Hk] = 0, for k < uq. For k — 2, the proof is the 

calculation given above in sections (b) to (g). For k = 3, the commutativity of H no+ i 

and H3 can be established using the Jacobi identity and the fact that [H no+ 2, H2} = 0. 

Similarly, one may successively show that [H no+ i, Hk >3 ] = using the Jacobi identity 

and the relations [H no+ £, H 2 ] = for t = 1, . . . , k — 1. The above method constitutes an 

alternative purely algebraic proof of the integrability of the XXX model. 



Notice that the formulae for 7 reflect our choice of normalization in (2.21); they need to be 
modified in another normalization, to take into account a rescaling of the Q n 's, corresponding to 
a rescaling of the spectral parameter. 
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5. The XYZ model revisited 

5.1. Proof of the general pattern for the charges 

We now present the proof that the charges of the XYZ model have the form ( |3.19|) . 
We will again proceed directly, calculating [B,f n (S)], for an arbitrary sequence S = C* 
with cluster C = {i±, . . . , i n } and pattern 7r = {ai, . . . , a n }. This commutator contains 
only terms of the form f n ±i(S'), where S' can be obtained from S by a few simple trans- 
formations, given below: 

<°_!S = {h ~ Ml, • . . , i n }{°o,a xa 1) ... ) a„} ) 

q a n + 1 _ f. ■ ■ j_ il{oi,...,0„+iXOn,O n +l} 

°°j n + l — I 4 l5 ■■•) ( n,%Tl/ , 

C^ — /,' ,', „ ,', , ,' \{oi)'"j°fe— 2j<*fc Xa fe _i,a fe + i,...,a„} 

<->^ a fc = \*1) • • • >*fc-l> 2fc+l 5 • • • i l nf 7 

« C — /,' /,', , ,*, , , ,' l{ a lv> a fc-i> a fc Xa fe+1 ,a fe+2 ,...,a„} 

^ a k>-> = \*1) • • • j *fc-l, *fc+l) • • • ) *nj , 

C a = /,', ,* . « (4, 4.\ i . « i \{ai,...,aj_i ) (a fe ,af) ! aj + i,...,a„} 



(5.1) 



+ <7. <7. 



The following conventions will be adopted: 

{oi,. . .,0, .. . ,a n } = 0, 
{ai, . . . , — afc, . . . , a n } = — {oi, . . . , a^, . . . , a n }, 

and 

/n(C°) =0, 

foe-") = - /an. 



(5.2) 



(5.3) 



The basic commutators of the form [jojOj+i, / n (<S)] may be evaluated similarly as in 
the XXX case. Summing up all these terms yields: 

[BJ n (S)] = b + (S)+b- + (S)+b + -(S)+b-(S), (5.4) 

where 

b a (S) = J2 6 e(S), a 6 {+,-+, +-,-}, (5.5) 

e=x,y,z 
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and 



b+(S) =-(»!- l)f n+ l(*l-lS) + infn+l(S<Tt n+1 ), 

*2 » n -l 

n-2 



^e (5) =Zl(A e ) / n+ l(«Sai eXO! )(1 — <^ 1)i2 _l) 

»1 »1+1 »1 

- (i n - l)(A e ) 2 / n+ i(5 cr a„^ (je o . e xa„)(l-5i n;in _ 1+ i) (5.6) 

n-2 

+ 2^^(A e ) /™-l( 5 (J °^ (T eX "j (7e )(1 - 5i j+1 ,i j +l) 



a . — ><r . a . 

j—2 l i l i l j + 1 

n-1 



+ ^2(h - 1 )( A e) 2 /n+l(5 (j a J .^ CTe CT eXa J )(l-^._ 1>ii _l), 

7=3 l J z j~ 1 % i 

K( s ) =2n(A e ) 2 /n-i(^r-5)5 il+ i, l2 - 2z n _i(A e ) 2 / n _i( l S— )S ln _ 1+lyin . 

With the exception of the third and fourth terms in 6" 1 (<S), all the terms in 
contain only expressions corresponding to ordered sequences. Fortunately, thanks to the 
factors arising from the coefficients g(S), the contributions due to disordered clusters will 
cancel" in summation over all possible clusters for a given pattern tt. Therefore, the 
commutator [B, F^ k ] can be expressed as a linear combination of the terms of the type 
f(S'), with S' being an ordered sequence. 

A simple inductive step finishes the proof. Assume that the n-th conserved charge is 
of the form ( |3.19| ). As we have just shown, the n + 1-th charge Q n +i, given by [B,Q n ], 
is a linear combination of the polynomials f(S). Now, because the commutator [B,Q n ] 
is translation invariant, it must in fact be a linear combination of elements of the spaces 
jF n _|_ 1 _ 2 fc ) ;. Therefore, the charge Q n +i is again necessarily of the form ( |3.19|) . 

5.2. Recursion relations for the F's 

Using the operations (5.6) (which naturally extend to any space J- n ,m)i one can derive 
recursion relations for the F's. It is convenient to represent this recursion as a transfor- 
mation of the triangle corresponding to Q n into a new triangle corresponding to Q n +i (cf. 
Fig. 2.). 



This can be seen by an argument similar to the one used in section (4.2.d) to demonstrate 
the cancellation of the r-terms in the XXX case. 
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p„+l 

1 n+1,0 




- n-3,2 



n-3,0 



Fig. 2. The action of the boost operator on the triangle representing Q n for the XYZ 
model. 



It is clear that the boost operator does not modify the structure of the triangle if n is 
even (i.e. the triangles corresponding to Q2m and Q2m+i differ only in labeling of points). 
For n odd, the action of the boost produces a new strip at the bottom of the triangle 
(therefore adding [{n + l)/2] new points). The left edge of the triangle representing Q n +i 
is given by: 



™+i 

± n+l-2fc,0<fc<[(n+l)/2] 



b- + {K_ 2{k _ 1) ,_ 1 )+b + {F 1 



n — 2k, k/^ 



(5.7) 



with 



F n+1,0 — u \ r n,0J 



The right edge of the Q n +i triangle is: 



(5.8) 



^+i- 2fc> o = &"(^n-2( fc -i),o) +b + (K-2k,o) +b + ~(F:. 2k , 1 ). 
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(5.9) 



The interior of the triangle can be obtained from: 

F n+l-2k,0<e<k = b ( F n-2(k-l),e) +b ( F n-2(k-l),£-l) 
+ b (F™-2k,l)+b (Fn-2k,l+l)- 

Finally, for n odd, the bottom edge of the triangle Q n +i is: 



(5.10) 



*&<«[(M-i)/2] = b-(F^) + b~ + (F^_ 1 ). (5.11) 

These relations provide a recursive way to calculate the conserved charges. As in the 
isotropic case, the basis {Q n } contains terms proportional to lower order charges; it is 
more convenient to express the results in another basis, denoted {H n }, obtained by taking 
appropriate linear combinations of the form ( p.22|) , in which lower order contributions are 
subtracted. 

5.3. Explicit form of the first few conserved charges 

We present below several of the lowest-order conservation laws of the XYZ model. 
The first non-trivial charge beyond the hamiltonian is: 

#3 = *3,0 = £}fo >< O-j+l) • *i+2- (5-12) 

The four-spin charge is: 

#4 = ^,0 + ^1 + ^0, (5-13) 

where 

jeA 

F 2,i=Z^ 2_, ( Aa ) g J g 3+2» (5.14) 

jeAaE{x,y,z} 

j'eAa6{i,i;,z} 

The five-spin charge corresponds to a triangle with the same structure as that of H4, 
namely: 

H 5 = F| i0 + F| 5l + F| )0 , (5.15) 
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with 

F 3,i=Z^ Z^ ( Aa ) e (^^+2^+3 + ^+1^+3), (5.16) 

j6Aa6{a;,i;,z} 

^5 \ ^ \ r / \ \2^abc~a~b ~c 

F 3,0 = -}^ 2_, ^ Xb > 6 a J a 3+l a J+2- 

jeAa^b^c 

The six-spin "monster" is: 

H 6 = F 6 6 + Fl A + F 2 6 2 + F 4 6 + Fl ± + F 2 6 , (5.17) 



with 



F lo =^2(((vj x CT i+i) x ^i+2) x dj- +3 ) x cr i+4 ) ■ or i+5 , 



jEA 

A f 



F 4,l = 2^ ZJ^ )(^ ^+2^+3^+4 - ^+2^+3^+4 



jeA a^b 



b~a ~b ~a ~a~b ~b 



1 ~b~a ~b ~a 

+ a,a, +1 a, +2 a,. 



a A CT- 1 i(T„- 1 oO" „• 



; w / l"j+2"j+4 " j " j+1" j+2" j+4 
;b 7^0, c-b ~a~b ~b 



I- A a A 6 (0 j cr, /+1 cr j+ 3Cr j+4 - <Tj<Tj+\<Tj+3<Tj+4,)\, 



FI2 =E E (^r^) 2 ^i + 3, (5-i8) 



ieAae{i,y,z} 



7716 V'Wx 2 \2\~a~6 ~6 -a 

^4,0 = 2^Z^ A & ~ A a) cr j cr i+l cr i+3 cr j+4: 

jEAae{x,y,z} 

F lo =E E ( A « + A ')( A ' + >Z)***j+i- 

The logarithmic derivatives of the transfer matrix can be expressed as linear combi- 
nations of the i? n 's: 

Q3 —H3, 

Q 4 =2H A - 2p 2 H 2 , 



Q 5 =6H 5 -4p 2 H 3: 

Qe =24H 6 - 32p 2 H 4 + 8p 4 H 2 . 



(5.19) 



with p 2 = X x + X y + X z . In the isotropic limit (X x — X y — X z — 1), the XYZ charges 
reduce to linear combinations of the XXX charges given by (|4.2| ) , which we denote below 
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as H xx to avoid misunderstanding: 



H 3 


=^3 


#4 


ttXXX j_ ttXXX 


H 5 


tjXXX ttXXX 
= H 5 ~ H 3 



(5.20) 

H 6 =H^ XX +4H? XX . 
From (5.19) we then recover the relations ( |4.66|) . 

6. The conserved charges in the XY model 

6.1. Generalities 

For completeness, we present in this section the results for the charges for the XY 
model: 

H XY = J>cr*crJ +1 + VKfi' (6- 1 ) 

As is well known, the XY spin chain is equivalent to a free fermion theory; the charges for 
this system have thus a particularly simple form, and they have been discussed by many 
authors [26-29]. Although the following results are not new, the derivation given here is 
different from the approach of [26-29] . An interesting feature of the XY case is that there 
exist two independent families of conservation laws. Both families persist when the model 
is perturbed by a perpendicular (i.e. in the z-direction) magnetic field. 

Since the XY model is just a special case of the XYZ hamiltonian, we can apply to 
it the apparatus developed in sections (3) and (5). The general formula ( |3.19|) greatly 
simplifies for the XY case. First, note that since A 2 = 0, the coefficient g(S) vanishes 
for any sequence containing holes. Furthermore, f n {<S) = for sequences with a± = z or 
a n = z or d2, • • ■ , a n -i = x, y. Thus, f n (<S) can be non-vanishing only for the sequences of 
the form: 

S nJ = i a f ' a j+l ' • • ■ ' a j+n-2 > a j+n-l } > ( 6 - 2 ) 

with a, (3 e {x,y}. Observe further that 

US^ 3 ) = e ef, (6.3) 



where 

^n,j = u j u j+l ■ ■ ■ u j+n-2 u j+n-l^ 
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et=a«a z 1+1 ...a z 1+ n-2°L n -n (6-4) 



and the coefficient e assumes values ±1 or 0. It vanishes when a 7^ /3 if n is even, or a = ft 
if n is odd. We also introduce the notation 



P af3 



£<3 ( 6 - 5 ) 



The charges ( ^.19| ) for the XY model take thus the simple form: 



[n/2]-l 

Q n = E ^-^o- (6-6) 

fc=0 

with F™_ 2k e J-'n-2k,o- This shows that the triangle corresponding to Q n collapses to its 
right edge. Note that the spaces J r n -2k,o are two-dimensional. A convenient basis is given 

by {e^_ 2 fc> e n-2fc} f° r n even, an d { e n-2fc' e n-2fc} f° r n °dd- As we wm shortly see, the 
family {Q n } contains only half of the integrals of motion of the XY model. 

6.2. The XX case 

Let us consider first the special case of the XX model with X x = X y . One finds two 
conserved quantities for each n, which define two families: 

#(+) = eT + el y n even, 



n ' n 



e x n y - el* n odd, 



and 

Tj(-) — P X V _ P V X 



"71 n 



e T + e l y n °dd. 



(6.7) 



(6.8) 



In particular, there are two two-spin hamiltonians; 

HP=52ofof +1 +*W +1 , (6.9) 

jeA 

^"^E^+i-^W+i- ( 6 - 10 ) 

jeA 

i?2 is J us t the XX hamiltonian, while H2 is a special case of the Dzyaloshinski-Moriya 
interaction [35]. 

Commutativity of H^ with both hamiltonians if^ can be easily verified directly 
by a short calculation. Furthermore, the boost operator, 



2%. % 
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has ladder properties: 



[B, fli+>] = (n - lX-ir'HH^ + ffW), (6.12) 



(valid for all n with the understanding that H± =0). Since the XX hamiltonian H 3 
is proportional to to the first logarithmic derivative of the transfer matrix Q2, it follows 
that H n ' can be expressed as a linear combination of the logarithmic derivatives of T. 
This shows the mutual commutativity of the Hn s. Note that the family {H n } can 
be obtained from the family {Q n } by a transformation of the form (|3.22| ). Clearly, the 
triangle corresponding to H n degenerates to a single point. In terms of the associated 
triangles, this transformation corresponds then to redefining to zero all the points on the 
right edge, except for the vertex F^ . 

It is clear from ( |6.12| ) that for even (odd) n, Q n can be expressed as a linear combi- 
nation of the Hm with even (odd) m < n: 

[n/2]-l 

Qn= E ^ )H n% (6-13) 

p=0 



The coefficients 7 satisfy the recurrence relation ( f4.65| ), where the coefficients (5 are now: 



Po = Pi = ( — 1)^ 1 ' ) ( n — !)• For example, up to additive constants, 



Q 3 


rr(+) 


Q 4 


= 2tf| +) -2^ +) , 


Q 5 


= 6^ +) +8^ +) , 


Qe 


= 24H^ +) + 40H { 4 +) + 16H^ +) 



(6.14) 



( |6.14| ) can be also obtained as a special case from the formulae ( p.!2[ )-( |5.17| ) for the general 
XYZ model. 

In contrast, the charges of the second family H n ' are not of the general form ( |3.19| ). 
Remarkably, the boost operator associated to if^ acts as a ladder operator for the second 
family too (as first noticed in [36]): 

[B, HM] = (n- l)(-iy(H n - + \ + H n ~_\), (6.15) 

with H[ ' = — 25^ eA cr|. Using (6.15), (|6.12j ), together with the Jacobi identity one can 
prove by an inductive argument similar to the one in section (4.2), that the charges of the 
second family commute among themselves and with those of the first family. 
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Note that under a spin rotation by tc/2 around the z-axis but only on even sites, i.e.: 

(6.16) 



a 2j -* a 2j , a 2j -► -a 2j 



a 2j + l ^ °2j + li a 2j+l ^ °2j+l 

the two-spin charge H 2 transforms into the staggered XX hamiltonian: 

jeA 

Both Hn and -f/A behave as scalars under a global spin rotation around the z-axis. 
Indeed, it is simple to check that all members of both families commute with the generator 
of such rotation, that is the z component of the total spin, 

S* = 5>J. (6.18) 

jeA 

Therefore, both families survive when the model is perturbed by a magnetic field term 
hS z . On the other hand, while H n is invariant under global parity transformations, i.e. 
a'j — ► — a?, H n changes sign and for this reason we qualify it as a pseudoscalar. 

6.3. The general XYh case 

We turn now to the analysis of the general XY model incorporating the magnetic field 
term. The defining hamiltonian is thus 

H XY + hS z = J2 A*<^; + i + VK+i + K- ( 6 - 19 ) 

jeA 

As in the XX case, we look for the charges in terms of the basis e^f . We denote by A n 
the span of {e n x , e n y } for n even, and the span of {e n y , e n x } for n odd (i.e. A n = J- n ,o)- 
Similarly, B n will denote the span of {e n y ,e n x } for n even, and that of {e n x ,e n y } for n 
odd. In addition, A\ — 0, B\ = span(5' z ). Let A = © n=1 A n and B = ® n=1 B n . We now 
introduce two derivations on A © B given by: 

6xY(P) = hH XY ,P], (6.20) 

5z(P) = ^[S\P], (6.21) 

for P e A © B. The action of these operators can be visualized on the diagram in Fig. 3. 
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Fig. 3. The action of the derivations Sxy , &z in the space A © B of the XY model. 

The derivation 5z acts "horizontally," i.e. 5z ■ A n — *■ <S n , 5z : £> n — ► v4 n) while 5xr 
acts "vertically," i.e. #xy : -An — ► «4 n -i © -4 n +i, and 5xy '■ B n ~ * &n-i © $n+i- The 
action of $xy on an arbitrary X E A® B can be written as: 



5 XY (X)=u(X)+d(X) 1 



(6.22) 



where tt : A n — »■ *4 n +i, w : £> n — > £> n+ i, and <i : *4 n — ► A n -i, d : B n 
calculation yields: 

MO=-A,(ef +1 + C +1 ), 

«(eD=Ax(er +1 +e^ +1 ), 
42 



B n -\. A simple 



(6.23) 



and 



d(e-)=-X x (el- 1 +e n y _ 1 ), 

d(e y n y )=X y (e y n x _ 1 +e x r l 1 ), (6.24) 

Similarly, 

mo =*(<r +«#'), 

*z(0=-M^ B + 0» (6-25) 

The kernel of u in the space ^4 n © i3 n is thus two dimensional (it is the span of 
e n V ~ e n X an d ^e^+Aye^), allowing for the existence of two families of integrals of motion. 
The explicit expressions for half of the XX charges are not affected by the anisotropic 



deformation of the hamiltonian: it is simple to check, using (|6.23| ), that 



H& = &-<%> (6.26) 

commutes with the hamiltonian, for any n. In order to find explicit expressions for the 
remaining charges, we write them in the form 

[n/2] l(n+l)/2] 

H i 2 > = X x e x n x + \ y e™ + J2 X n-2k + J2 Y n+i-2k, (6.27) 

fc=l fc=l 

with X n -2k e An-2k(B n -2k) and Y n+ i-2k e B n+ i-2k(A n +i-2k) if n is even (odd). The 

(2) 

condition [Hxy + hS z , Hn ] = gives then the following system of equations: 

u(Y n _ 1 )+5 z (\ x e x n x + \ y eiy)=0, 

(6.28) 
u(X n _ 2k ) + 5 z (Y n _ 2k+1 ) + d(X n _ 2k+2 ) =0, for k = 1, . . . , [(n + l)/2], 



which allow for a recursive determination of the X 1 Vs: using (|6.23|) , Y n -\ is fixed by the 



first equation; the second one, with k — 1, determines X n _2; then, the case k = 2 leads to 
y n _3, etc. Note that since the kernel of u in A m © £>m is two-dimensional, (6.28) does not 
fix the terms X n - 2k and Y n+ i- 2k uniquely, but only up to an admixture of lower order 
integrals of motion. In solving (6.28) one may thus impose additional constraints. One 
possibility is to demand that X m and Y m be orthogonal to the kernel of u in A m © B m . 
This yields 



n-l 



Hi 2) = KeT + \ y e\ y + J2 Ck(\ y e™ k - X x e y n y _ k ), (6.29) 
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where the coefficients C& satisfy 

1 

(X x + Xy ) 

with 



Ck = 7V~2 7~2^(H^x + Aj/)cfc_i - 2X x \yC k -2), (6.30) 



Cl = 7r~2 27 ' anc * c o = h 2 ; 2' (6.31) 



"•(Aaj — A v ) (A-y — A 3 

g 2T ' and C = T~^2 

(A X + Ay ) (A X + Ay J 

Alternatively, one may obtain a simpler solution by demanding that the number of terms 
X n _ 2 fc and Y n +\-2k in ( p.27| ) be minimal. This leads to following expressions, that have 
first been written in [26]: 

H^ = X x e x n x + \ v e y n y + r n _! + X n _ 2 , (6.32) 

where 

6.33 

-An-2 — ^xG-n-2 + A y e n-2' 

As in the XX case, it can be shown that the boost acts as a ladder operator. The 
recursive application of the boost, starting from the hamiltonian (6.19), generates a family 
of scalar conserved charges. For h = the n-th order charge of this family coincides with 
the n-th logarithmic derivative (2.18) of the XYZ transfer matrix, evaluated at X z = 0. 
The recursive application of the boost operator, starting from H^ produces a family of 
pseudoscalar conservation lawsJl3 

To end this section, let us mention that the degenerate case X y = of the XYh model 
(which could be called the Xh model) is known as the Ising chain. It is evident that this 
system inherits from the XYh model two distinct families of conserved charges. 

6.4- Higher order ladder operators 

An interesting feature of the XYh model is the existence of higher order ladder oper- 
ators [28, 36]. For simplicity we focus again on the XX case. The first moment of Hn 
acts then as a ladder operator whose application to a charge of order m generates a charge 
of order m + n — 1. For example, 



3 2i- 



jeA 



Note that for h = the charges of second family cannot be expressed in terms of the 
logarithmic derivatives of the XYZ transfer matrix, in the parametrization adopted here. Observe 
however that both families can be obtained from the XYh transfer matrix, written down in [37]. 
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acts in the following way: 

fo(+) rr(+)l _ tt(+) tt(+) 

[B 3 ,H 2 \-H 4 -H 2 , 

[Bf\Hi%] = (n-l)(H% 2 -H£\). { ] 

The first moments of the pseudoscalar charges transform scalar (pseudoscalar) charges into 
pseudoscalar (scalar) ones, e.g. 

acts as follows: 

[B { 2 -\H^] = (n- l)(If W - H™). (6.37) 

Since the commutator of two ladder operators has again the ladder property, any two such 
operators can be used to generate an infinite family of ladder operators. No such higher 
order ladder operators are known in the general nondegenerate case of the XYZ model. In 
particular, the first moment of the density of the three-spin XYZ charge does not have the 
ladder property. 



7. The su(M) invariant spin chain 

7.1. Formulation of the model 

In this section we consider the isotropic su(M) version of the XXX model in the 
fundamental representation. The hamiltonian of this system reads [38] : 

ffa = £W+i' (7-1) 

jeA 

where t?, a = 1, ...,M 2 — 1, are the su(M) generators in the fundamental representation, 
acting non-trivially only on the j — th factor of the Hilbert space &)jC N . For M = 2 
( [T.l|) is simply the hamiltonian of the spin-1/2 XXX model. It is convenient to choose the 
normalization so that t a are the su{M) Gell-Mann matrices: 

\t\t h ] =2if abc t c , 

(7.2) 

where f abc are the structure constants of su(M), and d abc is a completely symmetric tensor, 
non-trivial for all M > 2. 
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7.2. Expressions for the conserved charges 

Remarkably, the conserved charges for general M have the same structure as in the 
M = 2 case (the ordinary s = 1/2 XXX model). The expressions for the general case can 
be obtained from (|4.2|) by a simple substitution: 

_a ia ubc rabc (n q\ 

a -»• i e -»• / . (7.3) 



For example, 



Hs = £/ a6 % a *? + i^ +2 , (7-4) 



jeA 



h, = J2ir bp f pcd tp b J+1 t c J+2 t<* +3 + tjtj+i}- ( 7 - 5 ) 

jeA 
As for M = 2, we can define the spin polynomials f n (C), for an arbitrary cluster C = 
{ji,32,---,jn}, as follows: 

f tQ\ — fa,ia 2 b 2 fb 2 a 3 b 3 rb n - 3 a n - 2 b n -2 fb n - 2 a n - 1 a n j.a 1 j.a 2 jn» (7 g 1 } 

The general expression for the su(M) XXX conserved charges can then be written in a 
universal form, (corresponding to the Catalan tree in Fig. 1.), valid for all M: 

[n/2]-l fc 

H n = F n $ + 2_^ / J (^k,i^ p n-2k,e, (7.7) 

with 

F n , k = Yl f^- ( 7 - 8 ) 

The proof that the H n 's form a commuting family parallels the analogous proof in 
the su(2) case. First, we show that the identities (|4.8|) -( |i~T0| ) remain valid for su(M), but 
now with A = {t), . . . , tf 2 ' 1 } and B = {t l k , . . . , if 2_1 }, (with j ^ k) and with A x B 
understood as a vector with M 2 — 1 components given by 

(A x B) c = f abc A a B b . (7.9) 
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The demonstration of ( f4.8|) -( P~T0 ) is elementary and uses several identities for the tensors 
f abc and d abc M 

eabc cede rega _ _(M/2) f bd 9 

jabcjcde d ega _ _( M /2)d M9 , 

(7.10) 
rabc rabd -\/r scd 

ra.be rede reghihai radc rcbe reihihag 

The rest of the proof is exactly the same as in the su(2) case: from ( |4.8|) - (fO0 ) one derives 



( f4.16| ) and proves first [H2, H n ] = and then the mutual commutativity of all the H n 's. 

7.3. Relation to higher spin su(2) models 

The su(M) XXX system can be also regarded as a su(2) spin-s chain, with s = 
(M — l)/2. The hamiltonian (7.1) can be rewritten in terms of the variables S?, which are 



the spin-s su{2) generators, acting non-trivially only at site j: 

H = J2Gm(xj), (7.11) 

where Gm is a polynomial of order M 2 — 1, and Xj = SjSj +1 . In particular, the expressions 
for the first few polynomials Gm (obtained by expressing the matrices t a in terms of the 
standard spin-s matrices) read: 



Gs(x) =x + x 2 , 

G 4 (x) =x (11a; 2 + 4a; 3 ), 
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G 5 (x) =x + — (13a; -6a; -x), (7.12) 

G 6 (x) =x + — (14911a; 2 + 636a; 3 - 360x 4 - 32x 5 ), 
178435 

G 7 (x) =x + — — -(-6417a; 2 - 1713a; 3 - 8a; 4 + 19a; 5 + x 6 ). 
14094 



The proofs of the first three identities can be found in [39]. The fourth identity can be proved 
using the first two equalities as well as two additional well known identities: 

radc rcbe , rbdc race , rede rabc ^ 

/adc icbe , rbdc iace , rede iabe r. 

d +f d +f d =0. 
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It is convenient to use the basis provided by the projection operators Pjj, defined by 

pi . M = fi f^> < 7 - i3 > 

1=0 yi ^3 

where j — 0, . . . 2s, and 

y t = i[2l(l + l)- 8(8 + 1)]. (7.14) 

Acting on the states of the tensor product of two spins, Pf projects onto the states with 
total spin j. In terms of this basis, modulo constants, 

2s 

G M (X) = E(- 1 ) 7> (M-l)/ 2 (^)- ( 7 - 15 ) 

3=0 

Similarly, the density of a XXX conserved charge of order n in terms of the su(2) variables 
becomes a polynomial in n (matrix) variables Xj, . . . , Xj+ n -i. 

8. Equivalent representations of charges in spin chains 

8.1. su(2) spin chains in the Weyl representation 

It is well known that the su(2) generators can be represented as differential operators 
acting in the space of functions of two complex variables u, v satisfying |tt| 2 + \v\ 2 = 1: 

a z = ud u — vd v , 

a x = ud v +vd u , (8.1) 

a y = i{vd u — ud v ). 

This is often called the Weyl representation. The spin-s representation can be obtained 
by the restriction of the space of functions to polynomials of degree 2s. States of the 
spin-s chain with N sites can be therefore represented by polynomials of degree 2s in 2N 
variables Ui, v% (i — 1, . . . , N). The conserved charges H n become commuting differential 
operators of degree 2n in the variables Ui,Vi, d Ui , d Vi . For example, the XXZ hamiltonian 
in the Weyl representation is: 

H 2 =^2[(2 - K){ujV j+1 d Vj d Uj+1 +VjU j+1 d Uj d Vj+1 ) 

jeA (8.2) 

+ X z (ujU j+1 d U:J d Uj+1 + VjV j+1 d Vj d Vj+1 )]. 

48 



Higher charges can be obtained in the same manner, but the resulting expressions do not 
seem particularly transparent. 



Observe that in deriving the explicit form (|4.2| ) of the charges H n for the s = 1/2 XXX 
chain one makes use of both the su(2) commutator and the anticommutator a a a b + a b a a = 
5 ab . In the Weyl representation, this anticommutator is not encoded in (|8.1| ), but results 
from the restriction of the representation space to multilinear polynomials. It is therefore 



clear that the charges (|4.2|) cease to commute for spin chains with s > 1/2. Establishing 



the explicit form of the charges in this case is a much more difficult task. 

8.2. Representations of the XXX charges in terras of braids 

For both the s = 1/2 and the su(M) XXX models, the hamiltonian can be written, 
modulo constants, as: 

H = ^P jJ+1 , (8-3) 

where Pj y k are operators exchanging particles at sites j and k. The spin polynomials 
f n (C), where C = {i±, %2-, ■ ■ ■ , i n -i, in} can be rewritten in terms of the exchange operators 
as follows: 

This permits to express the higher-order hamiltonians for the XXX system as a sum of 
products of the exchange operators. 

Let hi = Pi,i+x be a nearest-neighbor exchange operator. It is simple to check that 
the ftj's satisfy the defining algebra of the braid group Bn, (where N = |A|): 

bibi+ibi = bi+ibibi+i, 

(8-5) 
bibj = bjbi, \i- j\> 1, 

and an additional constraint 

b\ = 1. (8.6) 

Note that any binary exchange (not necessarily nearest-neighbor) can be realized as a 
product of nearest-neighbor exchanges. Therefore, it follows that the spin polynomials 
( |S.-1| ) can be expressed in terms of the generators of B^. For example, if C contains only 
adjacent spins, 

f n {G) = [...[&!, 6 2 ],6 3 ],6 4 ],....,6 n ]. (8.7) 
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Hence the XXX charges can be interpreted in terms of equivalence classes (with respect 
to (8.6)) of braid polynomials (see Fig. 4. for an example). 





Fig. 4. A braid polynomial corresponding to the three-spin XXX charge 



9. Structure of conserved charges in the Hubbard model 

9.1. Introduction 

The Hubbard model (for which original references can be found in [40]) is usually 
defined in terms of fermionic operators, by the hamiltonian 

H=-2 J2 (4- a i+M + ii, s ^)+ 4 ^EKt- 1 /2)Kl- V2), (9.1) 



where 



n 3,s a j,s a 3\si 



(9.2) 



and U is some coupling constant. The sum over j runs from 1 to N, the number of sites. 
The operators a- s and a J;S are respectively the creation and the annihilation operator of 
an electron of spin s at site j, subject to the anti-commutation relation: 



t t X X 

Cdj s Ct>k,s' + dk,s'^j s = °j,kOs,s'i 



(9.3) 



and 



a ],8 a l,8> 



~ a k,s' a j,si 



O j j,sO j k,s f ^kjS' ™j,s* 



(9.4) 
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The hamiltonian can be rewritten in terms of spin variables by means of the standard 
Jordan- Wigner transformation: 



(9.5) 



a j,i = {° Z iV2-Vn){ T 1 T 2-Tj-i) Tj. 

{cTj} and {t"}, a G {x, y, x}, j G A, are two independent sets of Pauli matrices: 

K,r|]=0, (9.6) 



and the normalization used for a, is 



a± = i«±z^). (9.7) 



In this way, the hamiltonian H = H 2 is easily seen to be equivalent to (1.2). Notice that 



± 



for U = 0, the Hubbard model reduces to two uncoupled XX models and that for r = 



j 



and T? = h for all j, it reduces to the XX model in an external magnetic field. 

The first non-trivial conservation law was obtained by Shastry and found to be [6]: 

ieA (9.8) 

- U(af* v j+1 - ^oJ +1 )(r/ + r/ +1 )] + [*" r]. 

It is simple to verify that it indeed commutes with Hi- 

To conclude this short introduction of the Hubbard model, we point out that there 
exists another equivalent way to represent it: redefining the spin variables as a^ — ► S%j and 
Tj 1 — ► S%j+i, we obtain a lattice with 2A^ sites with next-to-nearest neighbor interactions 
and bond alternations: 

H 2 = J2 S+S~ +2 + S~S+ +2 + U[l + {-) j \S Z jS* +1 . (9.9) 

3 

9.2. The non-existence of a ladder operator 

We first demonstrate that for the Hubbard model there is no ladder operator B whose 
commutator with H 2 would produce if 3. The most general possible expression for B is 

[21]: 

B = J2^ b i+ c ^ ( 9 - 10 ) 

j 6 A 
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where bj and Cj are translation invariant densities, at most bilinear in (nearest-neighbor) 
spin variables. Enforcing the symmetry a <-> r and requiring B to reproduce the known 
ladder operator in the two limiting cases U — ► and r 2 — ► h gives 

b 3 =a*a* +1 + a]a y J+1 + t*t* +1 + rjrj +1 + |(a|r/ + ^ +1 r/ +1 ), 
Cj =0. 

The resulting S is thus the first moment of the hamiltonian, with the density of the latter 
symmetrized with respect to sites j and j + 1: 






This form of S is in fact generic to all the cases where it exists. 

With this candidate ladder operator, a simple calculation yields: 



[B, H 2 ] = -2z$>^ +1 aJ +2 - a>| +1 aJ +2 
jeA 



^K^i - ^K+i) W+i + r /)] + [°" ~ T l- 



(9.13) 



The sum has the same form as the above H$, except that the coefficient of the U term 
is here —1/2 instead of —1. This particular value of the relative coefficient between the 
U independent term and the linear one is of course crucial for the commutativity of H3 
with B.2- In other words, [H2, [B, H2]] does not vanish. This finishes the proof of the 
non-existence of a ladder operator for the Hubbard model. 

9.3. Higher order charges 

As already indicated in the introduction, higher order charges must be calculated by 
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brute force methods and they become rather complicated. For instance H4 reads:IIZJ 






'i+3 

jeA 



- £7(a 3 - ^ + i^ +2 + a 3 a ]+1 a 3+2 ){r J + r j+1 + r J+2 ) 

-zL(a x a y -n v n x Mt x t v -t v t x ) 
2 ^J a j' + l °3 a 3 + l^ T 3 T 3 + l T 3 T 3 + l> 

- Utfafa - aVa* +1 )(r* +1 Tj +2 - ry +1 r x +2 ) ( 9 " 14 ) 

+ U 2 (a]a^ +1 + vy j+1 ){T*T* +1 + 1)] 

+ [a <-»• r\. 
This charge has also been found in [41, 19]. Observe that these integrals of motion can 
be written much more compactly in terms of the tensor products of the densities of the 
XX charges. It is convenient to use a notation which differs slightly from the one used in 
section (6); namely, we define: 



l(+) 11 I yy rr x rr z rr z rr z /-r x _l_/-r^/-r 2 /-r 2 rr z rr y 

n n,j ~ e n,j T e n,j ~ a j a j + l a j+2--- a j+n-l cr j+n ^ a j a j+l' J j+2--- (J j+n-l <J j+m 

!,(-) — p x y _ P V X _^V rr z rr z rr y — a v n z rr z n z n x 

a n,j ~ e n,j e n,j — a j a j + l a j+2--- a j+n-l a j+n <7 j a j+l (J j+2--- a j+n-l a j+m 



(9.15) 



for n > 1, and 

h<$ = -a] , h[~) = 0, (9.16) 

and similarly 

9 { n ± J \r) = h^(a^r). (9.17) 

With this notation, the first three charges take the form: 

jeA 

H s =Tlhtt+9i; j ) + U{h£(g[y+g[%)+g£(h$ 
jeA 

H, =Y^[hi] + g$ + U {h$ (g$ + g[% + g[%), (9 . 18) 

jeA 

+ 9&\h® + h$ +1 + h[%) - h^ig^ + g£ + gtf +1 ) 



One could have hoped that the above operator B would still be of some use, at least for 
generating the correct terms in H4 (despite the fact that the coefficients cannot be expected to 
be right) in the commutator [B, H3]. But some terms in H4 are not produced in this way. 
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We have also obtained the explicit form of H$ which is 

H * =Efe ) + U {h{-J G,W + g[% + g[% + g{%) 

jeA 

+ h {+) (a { - ) +n { ~ ) + a { ~ ) + o (_) ) 

- «2j (S'lj-i + S'lj + Sij+i + 9i !j+2 )\ (9.19) 

+ t/ 2 {/i (_) fo (+) o (+) +0 (+) (+) +0 (+) (+) ) 
+ ^ 1^3,j K9i,j 9i,j+i -r ffij yi,i+2 + 9i,j+i9i,j+2) 

+ h (+) (a ( - ) a {+) + a (_) a (+) U 

+ n 2,j i^j-iyij+i + 9 2 ,j+i9i,j )} 
-U*{hbj ) (g$+gl%)}] + [h~g]. 

We now present a number of simple observations concerning the structure of the 
conserved charges, indicated by the above results. 

1- H n is composed of sum of terms of the form 

hUh%\..gWg%"\ (9.20) 

with ti > 1, i — 1, . . . ,p, where p is some integer and e, e', ... = ±. A parity of (9.20) is 
given by 

ee'...e"e'" = (-) n . (9.21) 

Furthermore, these terms must all be scalars, which means that in the product (9.20) 
pseudoscalar factors must occur in pairs .11^1 

2- With the normalization we have chosen for H 2, the relative coefficients of the various 
components of H n are ± times some powers of U. 

3- The term independent of U in H n is given by 

H n (U = 0) = J2lh ( ^+gi% (9-22) 

jeA 

with e = (— ) n . Terms depending upon U involve products of h and g factors of the form 

(9.20) with 

p 

^^=n-2r, (9.23) 



A pseudoscalar changes sign when the direction of the z-axis is reversed, h n is scalar 
(pseudoscalar) for n even (odd), while h n is scalar (pseudoscalar) for n odd (even). The same 
holds for g. 
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where r is a non-negative integer. These terms can be therefore grouped into classes, 
associated to the partitions of n — 2r, for r = 0, . . . , [n/2] — 1. 

4- The terms in H n associated to partitions of n (those for which r = 0) have a simple 
pattern. The number of parts p in the partition is related to the power of U in the 
coefficient of the corresponding term in H n . This power is simply p — 1. For example, the 
terms with r = which are linear in U are obtained from all possible decomposition of n 
into two parts of distinct 'colors' (the a and r factors), e.g. 

5 = 4 + (1) = 1 + (4) = 3 + (2) = 2 + (3), (9.24) 

(where terms referring to g factors are written in parentheses). Similarly the terms of 
order U 2 are associated with partitions of n into three parts of two colors. For instance, 
the terms appearing with a coefficient of absolute value U 2 in H$ are in correspondence 
with the partitions 

5 = 3 + (1 + 1) = 1 + 1 + (3) = 2 + (2 + 1) = 2 + 1 + (2). (9.25) 

But not all partitions into more than two parts are eligible (e.g., the partition 5 = (l)+2+2 
is absent). Some, but not all, are eliminated by the scalar and/or the parity conditions 
stated above. For instance, 5 = 1 + 1 + 1 + 1 + 1, with any arrangement of parentheses, is 
excluded by the parity and the scalar condition. 

5- The above correspondence is not one-to-one since more than one term is associated to 
a given partition. Consider a class of terms corresponding to a certain partition. Terms 
within such class are related to each other by a relative translation of their components. 
For the first two-component partition, i.e., n = n — 1 + (1), the different terms in the class 
correspond to the n — 1 positions of g[ k , k — j, ...,j + n — 1, relative to the /i n -i,j factor. 
For the other two-component partitions, there are also n — 1 terms in the class, with the 
difference that now the first term in the class n = n — k + (k) is gkj-k+i and the other 
terms are obtained by translating this g factor by one unit with respect to h n -kj, up to 
the site j + n — k. The different terms in classes labeled by partitions into more parts can 
be described similarly. 

It is more difficult to find a simple pattern for the terms associated to partitions of 
n — 2r, r > 1 because they are not universal, in the sense that they can be modified by 
the addition of a linear combination of the lower order charges if n _2 r - Nevertheless, the 
terms linear in U are described by the explicit formula presented below in section (9.5). 
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9.4- A- diagrammatic description of the conserved charges 

Formulated in words, some of the above rules may look a bit complicated. However 
they have a very simple diagrammatic description. To the densities h„ \ , g„ • , we will 
assign sequences of n symbols (corresponding to n successive sites of the chain starting at 
site j), as follows: 

ti n j — 0© • • • QjQjj 

ti n j = 00 • • • GG, 
g { n] =00 ■ ■ ■ 00, 

g { ~] =00 • • • 00. 
The product of two or more densities will be represented by writing the corresponding 
sequences of symbols, appropriately shifted, in successive rows, e.g: 

h (-)(+) _©0 

n 2,j 92,j+l ~ 00 • 

It is understood that the leftmost symbol in the top row is at site j. In this notation, 
(omiting the sums over j), the charges read: 

H 2 = 0© + 00 + U®, 



#3=00© + ^(| e + e |) + (^^ 

ffiffiffi ffiffiffi (TiiTiiTi 

if 4 =©©0© + U( ® 00 +®® +® ® ) 

' 

, © © © ) 

v 000 000 000 ' 

jji 00 00 00 \ 

^ 00 00 00^ 

u{ ® + ® + ® ) + u 2 ( ) + u\ © ) - u\®\ 

r, r,n 

00 

h 5 = eeeee + u ( 0eee + e0ee + ee0e + eee0 ) 



ffiffiffi (Ti(Ti(T1 ffiffiffi /7-\ /T\/T\ 

7-7- / tytfty tytyty tytfty w 

1 00 00 



jj i 00 00 00 00 



©e© eee eee 

r2 



+ t/ 2 ( +0 + ) + U 2 ( 00 + 00) 



- U "( ^ e + ®2 ) + (*~ T )- 



56 



9. 5. The explicit form of the term linear in U 

We now give the explicit expression of the term linear in U in the charge H n : 



tf» = £[<f +*£f 



jeA 

[n/2]-l n _2fc-l n-2 

/ \ ti — U fe 3 — I— nm I in — mn \ —I— I i I — I 

' ' " 7^ 

fc=0 m—\ e=o 

+ 0(U 2 ) 



+ U 2^ 1^ 1^^> h n-m-2k,j 9m,j-m-k+£+l\ (9.26) 



We now explain the origin of the three internal summations for the linear term. The first 
one, over k, keeps track of all the integers n — 2k whose two-term partitions are to be 
considered. The second summation takes care of all the possible ways the integer n — 2k 
can be separated into two parts. Finally, the sum over £ generates the different terms of 
a given class, that is the different translations of the g factor with respect to h, whose 
position is kept fixed. Although it may not be manifest, this expression is symmetric with 
respect to the interchange of h and g. 

To prove this result we have to show that, up to terms of order U 2 , (9.26) commutes 
with Hi- The proof is based on the following general commutation relations: 

[^°',*a ) ]=c*Si , .*a > ] 

ie\ (9.27) 

(=f) .*,(*) _i_J?f*,(T) _ »,(t) 



± 2i{fc&ij-i + h™ ld + S(fC! 1J+1 - h™ ld )}, 



and 



ieA (9.28) 

- ±z n n n,j9l,j n n,j9l,j+n-lf- 

In the first expression, we have introduced an operator S whose action on h n ( or g^ / is 
defined to be multiplication by 2 if n = 1 and 1 otherwise: 

6h?J=2h<$, Sh^=h^. (9.29) 

Exactly the same results applies for h — ► g. Notice that both type of commutators change 
the parity of the density on which it acts. 

These commutators are very easily demonstrated using the definitions (9.15), (9.16) 
and (9.17). 
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We now proceed to the evaluation of [H^ , Hn ]. Using (9.27), one gets 

[n/2]-l n -2k-l n-2 

jeA fc=o m=i e=o 



(9.30) 



x U — ) n + m + l \h { ~ > —h y ~> -\-r)h { ~> 

A U > l'n-m-2k+l,j-l " l n-m-2k+l,j ' un n-m-2k-l,j + l 

/ \n-\-m ( N7TI + 1 

un n-m-2k-l,j\ ym,j-m-k+e+l 

, /_Nm+l.(-J (-) _ (-) 

T v ^ " / n-m-2k,jiym+l,j-m-k+e y-m+l,j-m-k+£+l 

+ "fi'm-l,j-m-fc+€+2 _ "fim-lJ-m-fc+^+l-li- 

We first relabel the j summation index in each term in order to have all h factors at site j 
(e.g. j — 1 — ► j in the first term, etc.). Next, in the first, second, seventh and eighth terms, 
we relabel m as: m — 1 — > m. We then evaluate the summation over L All terms cancel 
two by two except for those at the boundaries of the £ interval. For the summation over 
??7,, a similar situation holds, except that now two values at both extremities of the interval 
contribute due to the 8 factors which are responsible for partial cancellations. Most of 
the resulting terms can again be canceled if in half of them we reshuffle the k index as: 
k + 1 — > k. Summing over k yields then 

which is exactly the same as — [if^ >-^n ]• This shows that, up to terms of order t/ 2 , 
[Hz, H n ] indeed vanishes. 



10. Concluding remarks 

The main concern of the present work is to exhibit the algebraic structure of conser- 
vation laws in integrable spin chains. For the XYZ model, this has led to a description in 
terms of simple multilinear polynomials in spin variables. The terms describing a charge 
H n can be grouped into classes corresponding to different types of clusters with prescribed 
number of spins and holes. The coefficients of the different terms satisfy recursion rela- 
tions. In the general anisotropic case, we have not been able to solve these relations. In 
two special cases however, for the XXX and XY models, explicit expressions for all the 
conserved charges can be found in closed form. The structure of charges that emerges for 
the Hubbard model is organized in terms of tensor products of densities of two XX models. 
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However, in the absence of a recursive structure we were able to prove this only for the 
zero-order and linear terms in the model's free parameter. 

For the XXX and XY models, the explicit construction of the charges provides an 
alternative proof of integrability, independent of the Bethe Ansatz or the transfer matrix 
formalism. Admittedly, this pedestrian proof is quite tedious and it does not even adress 
the issue of finding the eigenstates and eigenvalues of the hamiltonian. On the other hand, 
the proof for the Heisenberg model can be generalized in a completely straightforward way 
for the su(M) invariant XXX chain (which required the use of the nested Bethe Ansatz). 
One may therefore hope that there are other types of systems in which the unorthodox 
direct approach could be effectively used to prove integrability. 

One expects to find a similar pattern of conserved charges in integrable chains with 
non-trivial boundary conditions. In particular, for the XYZ model with nonzero bound- 
ary terms in the hamiltonian, the construction of higher order conserved charges can be 
achieved by modifying the densities of the XYZ higher order charges only near the bound- 
ary; in the "bulk," the densities of the conserved charges are not affected by boundary 
effects. A similar reasoning can be made for integrable chains with impurities. It is also of 
interest to find explicitly the deformation of higher order conserved charges in integrable 
models with a general quantum group symmetry. 

Another very interesting problem is to determine the structure of conservation laws in 
integrable chains with long-range interactions. The explicit expressions for the charges H3 
and H4 in the Haldane-Shastry model [42] are analogous to the corresponding formulae for 
the Heisenberg chain. Indeed, one can regard the Heisenberg model as a limiting case of a 
long-range model. It is thus natural to expect the structure of the charges in the models 
of the Haldane-Shastry type to be similar to the XXX case. 

In a different perspective, the structure of the higher conserved charges for the models 
analyzed here may provide new insights into the problem of testing integrability for general 
spin chains. The integrability of a quantum chain is usually demonstrated rather indirectly 
by showing that the model can be solved by the coordinate Bethe ansatz, or that the 
hamiltonian can be derived from a commuting family of transfer matrices related to the 
Yang-Baxter equation. But these are only sufficient conditions for integrability. Moreover, 
testing these sufficient conditions is often not easy and systematic. Alternatively, one may 
look directly at the existence of higher order charges, with a general structure similar as 
in the basic models considered in this work. A heuristic integrability test based on this 
approach will be considered in a forthcoming publication. 

59 



For higher spin chains, anisotropic su(M) chains, or models with underlying algebra 
different from su(M), the determination of the explicit form of the conservation laws 
remains a difficult challenge. As we have stressed in the introduction, the problem is 
not only with the actual computation - here one may be helped, to a certain extent, by 
effective computer algebra programs - but in finding a pattern in the huge amounts of data 
that emerge out of these calculations. In our view, integrability reflects itself into such 
structural patterns. 

Appendix A. Mastersymmetries and hamiltonian structures in classical soliton 
theory 

A.l. Symmetries and mastersymmetries 

We consider continuous systems subject to the evolution equation 

<p t = K 1 (<p), (A.l) 

where ip collectively describes the independent fields under consideration. K\ stands for 
a vector field whose evaluation at the point ip yields the scalar quantity Ki(ip). The Lie 
derivative in the direction of Ki, noted L^ 1 , is the same as the time derivative along the 
integral curves of the above evolution equation. Thus, the time derivative of a generic 
tensor field A, that can depend explicitly upon time, is simply dtA + Lk x A. A is said to 
be invariant with respect to (A.l) if its time derivative vanishes, 

d t A + L Kl A = 0. (A.2) 

A vector field satisfying (A.2) is called a symmetry of (A.l). The Lie derivative of a vector 
field is defined as 

L K A = [K, A] L = A'[K] - K'[A], (A.3) 

where If' (<£>)[>!] denotes the Frechet derivative of K at the point ip in the direction of A: 

K'[A] = ^-K^ + eA)\ e=0 . (A.4) 

If (A.l) is one member of an integrable hierarchy, then there exists an infinite number of 
time independent vector fields K m all commuting together: 

[K n ,K m ] L = 0. (A.5) 
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All K m> i are symmetries of (A.l). A mastersymmetry is defined to be a vector field r 
whose commutator with any K n lies in the commutant of K n [43, 44, 45] : 

[[r,K n ],K n ] L = 0, Vn. (A.6) 

A. 2. Hamiltonian structures in integrable systems 

The evolution equation (A.l) is said to be hamiltonian [46,47] if it can be written in 
the form 

dip 
where the Poisson bracket is defined in terms of the differential operator P by 

{<p{x)Mv)} = P(x)8(x-y) (A.8) 

and 

A (matrix) differential operator P defining a Poisson bracket, satisfying the usual require- 
ments of antisymmetry and the Jacobi identity, is said to be hamiltonian. Take for instance 
the KdV equation 

ut = u xxx + 6uu x . (A. 10) 



It can be written as a hamiltonian system in two different ways [48] 

5 JL 3 - a S /V„.3 l ..2 
Su 



u t = P l"-T 1 = d *fa f( u3 ~ \ u D dx 



(A.ll) 
P 2 I^11 = UQi + 4udx + 2u x )^- I u z dx. 



S 4 1 = m + ^d x + 2u x )f I 2 
ou ou 

We stress that Pi and Pi define two distinct Poisson brackets (cf. (1.10)). (Notice also 

that here the subscript for the conserved charges stands for half the usual degree of the 

density: deg u = 2 deg d x =2.) Any equation in the KdV hierarchy can be written in the 

form 

SH n +i „ SH n , 

u t = Pi — - s =P2^—- ( A - 12 ) 

ou ou 

This translates into the Lenard scheme (1.11) for calculating the conserved charges. Such 
a system is said to be bi-hamiltonianllB and it is characterized to a large extent by its 
recursion operator 

R = P 2 P^ (A.13) 



More precisely, it is also required that any linear combination of Pi and P2 must be 
hamiltonian. 
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(the adjoint of the Lenard recursion operator). It is not difficult to show that if both Pi 
and Pi are hamiltonian, then the infinite sequence of operators 

P m = R m ~ 2 P 2 (A.14) 

are also hamiltonian. 

A. 3. Hamiltonian master symmetries for the KdV and NLS equations 

A mastersymmetry not only maps (via L T ) symmetries into new symmetries, but also 
any invariants (conservation laws, hamiltonian operators, mastersymmetries) into other 
invariants of the same type. In particular, there is an infinite number of mastersymmetries. 

An almost universal symmetry in integrable systems is the scaling symmetry and this 
can be used as a convenient structural organizing tool for the description of mastersym- 
metries [44, 45]. Let tq be the time-independent part of the vector field associated with 
this symmetry (an example is given below). The infinite set of mastersymmetries can be 
denoted by {r n } where in terms of the scaling, r n has degree n. The mastersymmetries 
satisfy the algebra 

[r n , Trn\L = c(m - n)r n+m , (A. 15) 

where n, m > and c is some constant. We also have 

[T n ,K m ] L = a n , m K n+rn (A. 16) 

(with a n}7n = constant). These relations show that the whole structure of the hierarchy can 
be extracted out of a degree 1 mastersymmetry r±. Actually, it is natural to expect that 
for an evolution equation of the form (A.l), the existence of a degree 1 mastersymmetry 
is sufficient to ensure the existence of an infinite number of mastersymmetries, a definite 
characteristic of an integrable system£2l In terms of n, the conserved charges can be 
generated recursively by 

L Tl H m = H' m [r x ] = H m+1 . (A.17) 



This point of view has been mainly advocated by Fuchssteiner. A related conjecture is 
that integrability is ensured by the existence of one non-Lie point symmetry [45] (section 4 of 
this reference explains the connection between these concepts). For integrability tests based on 
symmetries, see [49]. 
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A mastersymmetry is called hamiltonian if it can be written in the form 

r = P^ (A.18) 

for some hamiltonian operator P and some integral B. For bi-hamiltonian systems with 
scaling symmetry, all mastersymmetries appear to be hamiltonian with respect to an ap- 
propriate P m . 

Again, the KdV equation offers the simplest illustration of this feature. We first note 
that this equation is invariant under the following scaling 

u(x, t) -»■ e 2 u(ex, e 3 t). (A.19) 

This symmetry is generated by the time-dependent vector field 

r = 2u + xu x + 3tu t: (A.20) 

whose time-dependent part is the mastersymmetry To- It is not difficult to check that 

L T0 R = 2R. (A.21) 

An infinite number of mastersymmetries can be generated out of tq by the recursion 
operator R [44]: 

Tn = R n To = (02 + Au + 2u x d- 1 ) n (2u + xu x ). (A.22) 

These r n satisfy the algebra (A. 15) with c = 2. In particular 

T\ = x(u xxx + 6uu x ) + 8u 2 + Au xx + 2u x (d~ 1 u). (A. 23) 

Then, denoting 

B = \ I xudx, (A.24) 

all r n are seen to be hamiltonian, each with respect to a different hamiltonian operator 

[50]: 

r-i = Pi-f-, 

ou 

r, SB 
TO =^2^, 

OU 



ou 



-p 5B 

OU 
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SB (A.25) 



Since 

r n — — dx = H n+m , (A. 26) 

the above relations can be reexpressed in the form 

{B, H n } m = H n+m _ 2 (A. 27) 

(where { , } m is defined in terms of P m ). This shows that a KdV raising ladder operator 
(mapping H n to H n+ i through some Poisson bracket) is necessarily defined in terms of a 
composite hamiltonian operator, namely 

^3 = P 2 Pr 1 P 2 . (A.28) 

Equivalently, T\ is hamiltonian, but with respect to the composite operator P 3 . With 
respect to P 2 , B is not a ladder operator. On the other hand, the first moment of the 
hamiltonian H 2, that is ^ J xu 2 dx, cannot serve as a ladder operator with respect to P 2 
because 

P2j~\ I xu 2 dx = P 2 (xu) (A.29) 

is not a mastersymmetry. 

For the nonlinear Schrodinger equation, the analogue of (A. 27) also holds, where now 



B= xcf)*(f) dx, (A.30) 

and the Poisson bracket defined by Pi is 

{r(x),<P(y)} = 5(x-y), (A.31) 

whose quantization yields the usual form of the quantum Schrodinger equation. This B 
acts as a lowering ladder operator: 

{B,H n } 1 = H n _ 1 . (A.32) 

For the quantum version of these theories, the implication of these results is clear. 
The quantum nonlinear Schrodinger equation is obtained from the quantization of the first 
hamiltonian structure while the quantum KdV equation is obtained by the quantization of 
the classical KdV second hamiltonian structure. As already stressed in the introduction, 
for quantum systems the bi-hamiltonian property is lost. The fact that for both cases no 
classical ladder operator can be defined in terms of the hamiltonian structure appropriate 
to its quantum version, strongly suggests that no quantum ladder operators exists for 
the corresponding quantum theories. This conclusion is supported by a direct analysis 
presented in appendix B. 
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A. 4- Hamiltonian master symmetry of the Landau- Lifshitz equation 

The situation described above for the KdV case applies to most bi-hamiltonian sys- 
tems. However, there exist exceptional systems for which the mastersymmetry of degree 
1 is hamiltonian with respect to one of the two basic hamiltonian operators -Pi, 2- This is 
exactly what happens for the Landau-Lifshitz equation, to which we now turn. 

The canonical equation of motion defined by the XYZ hamiltonian 

jeA 



is 

da? 



-j*- = [cr%, H] = 2i\ h e abc a c j {a b j _ 1 + a b +1 ). (A.34) 

Its continuous classical limit is obtained from the substitutions [51]: 

<r?^S a (x), a? +1 ^S a (x + A), 

A 2 (A.35) 

X a ^l + —J a , t^-iA~ 2 t, 

with A being the lattice spacing, and it reads 

dS a 



dt 



2e abc S c (J b S b + S b xx ). (A.36) 



(A. 36) is called the Landau-Lifshitz equation. It can be written in hamiltonian form by 
means of the following Poisson bracket 

{S a (x), S b (y)} = 2e abc S c (x)5(x - y), (A.37) 

and the hamiltonian 

H=\ f{SlS a x + J a S a S a ) dx. (A.38) 

It follows from the field equation that S a S a is time independent; we can thus set 

S a S a = 1. (A.39) 

To proceed, we introduce a more compact notation: 

S-S = S a S a , J = diag{j\j 2 ,J 3 } , (SA ) ab = 2e abc S c , (A.40) 
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in terms of which the equation reads 

^£=SA (S xx + JS) = SA 5 -§. (A.41) 

at ob 

The mastersymmetry of degree 1 for this system is [24] 

n = x(S AS XX + S AJS) + S AS X . (A.42) 

Quite remarkably, it is hamiltonian with respect to P± = SA: 

ri = 5A^, (A.43) 

with B given by 

B = J x(-S x -S X + S-JS) dx, (A.44) 

which is exactly the first moment of the hamiltonian. 

The Landau-Lifshitz equation is integrable and it has an infinite number of conserved 
integrals H n (the density of H n contains n S- factors) . The corresponding symmetries are 
defined by 

Kn = SA ~w- (A - 45) 

These vector fields can be generated recursively from T\ by 

[ri,K n ] L = const K n+1 , (A.46) 

which is equivalent to 

[SA 6 -^,SA 6 -^] L = const S A ^±i. (A.47) 

ob ob ob 

Using the standard isomorphism between the Lie commutator and the Poisson bracket 

[s A % ' 5 A lt h = SA Js {{B ' Hn}) ' (A - 48) 

we end up with [24]: 

{B, H n } = const H n+1 . (A.49) 

Therefore, the quantum ladder operator of the XYZ chain model survives in the continuous 
classical limit. 
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A. 5. Higher order master symmetries of the classical XX model 

Although higher order ladder operators have not been found for the quantum XYZ 
model, they have simple form in the degenerate case of the XY model, as we have seen 
in section (6). One can therefore expect the existence of an infinite family of hamiltonian 
mastersymmetries {r n } in the classical continuous limit of the XY model. This is indeed 
the case, as we now show. To simplify discussion we focus on the XX case. 

Using the Jordan- Wigner transformation (9.5), the XX conservation laws (|6.7| ), ( |6.8| ) 
can be written, modulo constants, in the form: 

#n ±} = J>/+n-i°J T a+aj+n-!, (A. 50) 

jeA 

where the fermionic variables satisfy the anticommutation relations (9.3)- (9.4) (leaving out 
unnecessary subscripts) . To obtain the classical continuum limit we represent the fermionic 
variables in terms of Grassman (anticommuting) functions as follows: 



a 



+ — £ + 
j+n 



£, + (x + n5), a j+n = £ t (x + n5). (A. 51) 



The Poisson bracket is fixed to be: 

{£(*), t(y)} = {t(x),av)} = S(x - y), (A.52) 

ti(x),Z(y)} = {t(x),t(y)} = o- (A.53) 

The corresponding classical conservation laws become then (modulo a linear combination 
of lower order charges) : 

H& = |[(dr ^^ T ^(d^Oldx. (A.54) 

Integrating by parts, one obtains 

H&> = |[(1 ± (-V^dr'n^dx, (A.55) 

from which one can see that -f^n+i an d H^ n are zero. In other words, even though the XX 
model has twice as many charges as the nondegenerate XYZ model, half of these charges 
disappear in the continuous classical limit. 
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The hamiltonian becomes 

H 2 = \Hi +) = f^ x dx, (A.56) 

giving the canonical equations of motion: 

& = &, & = C (A-57) 

The degree zero mastersymmetry is easily found to be 



xii + c 



ro=(^ + i% + )- (A-58) 



Indeed, its Lie commutation with the vector field 



K = [ r+ ) (A.59) 



gives [to, K]l = 0. Higher order mastersymmetries are easily found as 

T "=( a »C/+))' (A - 60) 

They satisfy the algebra 

[r n , r m ] = (m - n)T n+m . (A. 61) 

Consider the hamiltonian structure P\: 

Pl= (i o)' (A - 62) 

whose components are defined by 

P,6(z-v)-( UW'tW itW'^iv)} \ (A 63) 

All the mastersymmetries in (A. 60) are hamiltonian with respect to Pi, with the quantities 

B n+2 = fxt+dZ+^dx, (A.64) 

playing the role of the hamiltonians: 

r n = Pi ( ? ) B n+2 . (A.65) 
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A. 6. Remark on the classical Heisenberg chain 

We conclude this appendix with a remark on the classical lattice Heisenberg model. 
It is interesting to note that the quantum ladder operator does not survive when only 
the classical limit is considered. Consider for simplicity the isotropic XXX model, whose 
classical version is defined by the hamiltonian [52, 53, 54] 

# = J>(1 + S^ +1 ), (A.66) 

j'eA 

with the constraint 

SfSf = 1. (A.67) 

This expression for the hamiltonian may look somewhat surprising at first sight. The naive 
choice would have been ^Zj € j^SfSf +1 , but it turns out to be non-integrable. In particular, 
the candidate H3 having the structure of its quantum relative does not commute with this 
naive form of the hamiltonian. Furthermore, the canonical equation of motion it induces 
has no zero-curvature representation [47]. On the other hand, recall that the classical 
limit refers to large values of the spin, and that the integrable higher spin version of 
the quantum Heisenberg model involves higher powers of SfSf +1 , i.e. for spin s, it is a 
particular polynomial of degree 2s in SfSf +1 [55]. It thus appears that in the classical 
s — ► 00 limit, the logarithm captures the essential integrability aspect of this polynomial. 
The lattice Heisenberg hamiltonian (A.66) with the Poisson bracket 

{S%,S]} = 2e abc 8 ij S c j (A.68) 

defines the classical version of (A. 34). Somewhat surprisingly, there is no ladder operator 
defined with respect to this Poisson structure. To demonstrate this, we introduce 

jeA 

where / is a function to be determined. By requiring {B, H} (with H given by (A.66)) to 
be independent of j, f(X) is forced to be ln(X). The candidate H% so produced is of the 

form 

, qa qb qc 

H 3 = V e °tea J o J - +1 a J - +2 

jeA ( 1 + S j S j+i)( 1 + ^+1^1+2) 
whose Poisson bracket with H does not vanish. 
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Appendix B. No-go theorem for the existence of a ladder operator for contin- 
uous integrable systems related to the XYZ model 

B.l. The nonlinear Schrodinger equation 

The first four conservation laws of the quantum nonlinear Schrodinger equation are 

H = fq+q dx, 

H 3 = [ {y+ x V x + 3V+y + W)dx, 



(b.i; 



and the defining commutation relations are 

[*+(x), *(y)] = S(x - y) , [*(x), *(y)] = 0. (B.2) 

For definiteness, we treat the repulsive case, and the coupling constant has been rescaled 
to one. 

The problem that we now consider is to try to find an operator B such that, modulo 
lower order charges, 

[B,H n ] = c n H n+1 for n=l,2,..., (B.3) 

where the coefficients c n 's are not zero for all n greater then some uq. Restricting our 
attention to the four first integrals of motion, which contain terms with at most two 
creation and two annihilation operators, allows us to reformulate the problem in the two- 
particle sector without being restrictive. 

Recall that a complete set of eigenfunctions in the iV-particle sector is given by (see 
e.g. [56]) 

|Ai,...Aat) = -= / X n{xi,...,x n \X 1 ,...,X n )^ + {x 1 )...^ + (x n )\0), (B.4) 

where 

EN 
., ^, , Jjr .,, ..,-, .^ v .j •.,,,. -=i*" A ^. (B.5) 

P j>k 
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P is a permutation of the set (1,2,..., N). In the two-particle sector, the expression of the 
conservation laws -Hi, 2,3 is (up to irrelevant multiplicative factors) 

H 1 =d 1 + a 2 , 

H 2 =d$ + d$-26(x 1 -x 2 ), (B.6) 

£3 =df + di- 3S(x! - x 2 )(d 1 + d 2 ), 

with di = d Xi . 

We now look for a differential operator b such that 

[H u b] = Cl H 2 , [H 2 ,b] = c 2 H 3 , (B.7) 

for some constants c 2 and c\ 7^ 0. It is simple to see that the first relation is satisfied by 
the following choice for b: 

x\d\ + x 2 d 2 - (xi - x 2 )S(xi + x 2 ) (B.8) 

(this gives c\ = 1). Actually this first commutation relation determines b only up to a 
term of the form 

f(x 1 -x 2 )M(d 1 ,d 2 ), (B.9) 

where / is an arbitrary function of the difference x± — x 2 and M is a general polynomial 
in the derivatives, which does not contain any Xi dependent terms. The commutator of 

b = x x d\ + x 2 d\ - On - x 2 )8( Xl + x 2 ) + f{ Xl - x 2 )M(d u d 2 ) (B.10) 

with H 2 yields 

[H 2 ,b] = 2(df + dl)M-45(d 1 +d 2 )+2(djf)M+2(d 1 f)M(d 1 -d 2 ) + 2f(M5-5M), (B.ll) 

where 5 is a shorthand for S(xi + x 2 ). We want to write the rhs of the last expression in 
the form 

c 2 H 3 = c 2 [df + di- SS(x! - x 2 ){di + d 2 )]. (B.12) 

This leads to the following equation for the functions / and M: 

2(8ff)M + 2(d 1 f)M(8 1 - d 2 ) + 2f(MS - SM) = (c 2 - 2)(df + 8%) 

(B.13) 
+ (4-3c 2 )S(x 1 -x 2 )(d 1 + d 2 ). 
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First suppose that c 2 7^ 2. Given that the higher order operators on each side of this 
equation must be the same, under the assumption that (dif) 7^ 0, it forces 

2(d 1 f)M(d 1 - d 2 ) = (c 2 - 2)(d 1 3 + dl). (B.14) 

Since there is no xi dependence on the rhs, (dif) has to be a constant. With 

/ = ai(xi - x 2 ) + a 2 , (B.15) 

where by hypothesis ai 7^ 0, M is found to be 

M = ^-(d* + d 3 2 )(di - d 2 )-\ (B.16) 

01 — 1 

Equation (B.13) reduces to 

2f(MS - 5M) = (4 - 3c 2 )5(di + d 2 ). (B.17) 

With the expressions just obtained for / and M, the last equation is not satisfied. There 
is thus no solution when c 2 7^ 2 and (dif) 7^ 0. A similar argument shows with the 
assumption (dif) = there is no solution either. With c 2 = 2, (B.13) becomes 

2(dff)M + 2(dif)M(di - d 2 ) + 2f(MS - 5M) = -25(x t - x 2 )(d x + d 2 ), (B.18) 

for which it is easily seen that again there is no solution. This demonstrates the non- 
existence of an operator b that could satisfy (B.7). In other words, (B.3) is impossible for 
Ci 7^ 0. A similar argument, making use of the Jacobi identity, excludes the possibility 
that ci = . . . = c no _i = with c no 7^ in (B.3). 

B.2. The quantum KdV equation 

In terms of the quantum extension of the classical KdV field, denoted by T, the 
quantum KdV equation reads [23]: 

d t T = [T,H] , H=^-idw (TT)(w), (B.19) 

2m J 

The parentheses stand for the usual normal ordering [57] i.e. 

1 f dx 
(AB)(w) = & ——A(x)B(w), (B.20) 

2-ki /,„ x — w 



IV 
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and the subscript w indicates that the contour integration circulates once around the 
singular point w. The second Poisson structure is now replaced by its quantized form, 
which we take as the radial operator product expansion [58] (the global form of the Virasoro 
algebra [59]) 

TV Y7Y ^ c/2 2T(w) 2d w T(w) 

T(z)T(w) = - - I + - r^ + + ..., (B.21) 

(z — ty) 4 (z — w) z z — w 

where the dots indicate irrelevant regular terms and c is the central charge, a free parameter 
of the quantum theory that cannot be scaled out. The commutator in (B.19) translates 
into an operator product expansion 

[T(z),H] = ^-idw T(z)(TT)(w). (B.22) 



This equation can be written as 



(! " c) & 



d t T = v J d 6 z T - 3d z (TT). (B.23) 



The first three conservation laws are [23] 

H 1 = <bT dz, 



H 2 = <f(TT) dz, (B.24) 

H 3 =j[{T{TT)) - Q±A{p M Td z T)\ dz. 

Following the same strategy as in the nonlinear Schrodinger case, we look for an 
integral B satisfying (B.3). The first commutation relation is satisfied with the choice 

B = J> z(TT) dz + S , (B.25) 

where Bq is in the commutant of Hi. With Bq = 0, it is simple to check that [B, H2] is 
not proportional to H%. 

Can we find a non-zero So that commutes with Hi? Clearly Bq = § z n A dz is not in 
the commutant of Hi for any n > 0. If n = 0, the only dimensionally possible choice for 
a local expression for A is proportional to d z T, which makes So vanish. The remaining 
possibility is to have a non-local expression for A. Because the pole term in the operator 
product T(z)A(w) has residue d w A(w), we see that this will indeed commute with Hi. 
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However, in the commutation with H2, it will produce non-local terms; since there are no 
such terms in if 3, this last possibility is ruled out. 

These arguments can be made even more explicit at c = —2 [60]. For this value of the 
central charge, T can be represented by the bilinear 

T=(H>), (B.26) 

where <p and ip are both fermions of spin 1 with operator product expansion 



T2 , 1p(z)(j)(w) = 

(z — w) z (z — wy 



</>(z)i/>(w) = -_ , V(^)0H = 7 w- (B.27) 



In terms of these fields, the infinite set of conserved quantities is 

H k/2+1 = (f dz(^ k ^)(z), (B.28) 



where 4>^ k > = d k (j). For k even these integrals can be reexpressed in terms of T as follows 
[60, 61]: 

on— 1 r < — 

H n = ——j)dz(T n )(z), (B.29) 

where 

{f n ) = (... (((TT)T)T) ...T) (n factors). (B.30) 

From 

(Tr)(z) = §(0'ty + W'')(*), (B.31) 

where a prime stands for a derivative w.r.t. z, we have 



B = B l + B with B 1 = (b zcj)"ipdz. (B.32) 

A simple calculation gives 

[# 2 ,Bi] = (2ttz) iw[^"^'" + (5) V>]M dw, (B.33) 

which is clearly not proportional to H3 = § <j>^ip dw because 

[0V" + (B V]#0*(-)- (B.34) 
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At c = —2, all conservation laws are bilinear in the two fermionic fields so Bq is also neces- 
sarily bilinear (otherwise non-bilinear terms would be generated from single contractions) . 
The most general non-local and bilinear form for Bq is 

B = a I z(d- m (j))(d™ +2 (j)) , a = const. (B.35) 

This is indeed in the commutant of H\. By integration by parts, this can be rewritten in 
the form 

B = a(-) m I z<j>"i/> dz + a(-) m (m + 2) 1 0ty dz. (B.36) 

The first term is proportional to B\, while the second one is simply .f/3/2, which commutes 
with all H n . Thus, with the most general form for B, we must still conclude that [H2, B] 
cannot be proportional to if 3. 
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